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O ! ABSTRACT 

O 

We present a study of the trends in luminosity, linear size, spectral index, and redshift 
of classical double radio sources from three complete samples selected at successively 
fainter low radio-frequency flux-limits. We have been able to decouple the effects of 
the tight correlation between redshift and luminosity (inherent in any single flux-limited 
f-^ \ sample) which have hitherto hindered interpretation of the relationships between these 

four source properties. The major trends found are that (%) spectral indices increase with 
linear size, (ii) rest-frame spectral indices have a stronger dependence on luminosity than 
on redshift except at high (GHz) frequencies, and that (Hi) the linear sizes are smaller 
at higher redshifts. We reproduce the observed dependences in a model for radio sources 
(born throughout cosmic time according to a radio-source birth function) whose lobes 
are fed with a synchrotron-emitting population from compact hotspots, and which suffer 
inverse Compton, synchrotron and adiabatic expansion losses. The magnetic energy 
density within each hotspot is proportional to the jet-power, and synchrotron losses 
suffered in the hotspot mean that the energy spectrum of the emitting particles fed to 
the lobes is governed by the jet-power. The axial ratios of radio sources in our model 
increase as the sources age, and axial ratios are higher in sources with higher jet-power. 
In simulating the basic observed dependences, we find that there is no need to invoke any 
systematic change in the environments of these objects with redshift if the consequences 
of imposing a survey flux-limit on our simulated datasets are properly included in the 
model. It is also necessary to include appropriate energy loss mechanisms (such as 
the effects of the cosmic microwave background and feeding the lobes from a compact 
hotspot) which cause decreasing luminosity through the life of a source. Although our 
study has broken the luminosity-redshift degeneracy, we present evidence that for such 
studies there is an unavoidable "youth-redshift degeneracy" , even though radio sources 
are short-lived relative to the age of the Universe; it is imperative to take this into 
account in studies which seemingly reveal correlations of source properties with redshift 
such as the "alignment effect" . 

Subject headings: radio continuum: galaxies — galaxies: evolution — galaxies: jets - 
galaxies: active — quasars: general 
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1. Introduction and scope 

In order to explore the physical mechanisms by 
which classical double radio sources are governed, ob- 
servations of samples of such objects are needed to dis- 
cern the salient relationships between different source 
properties, such as spectral index or luminosity. 

Even observations of a complete sample will in- 
evitably be convolved with the cosmology of the Uni- 
verse as well as with the consequences of the sample 
selection process. We first describe (in §^|) how to 
avoid compounding the latter problem and giving rise 
to spurious, or hard to interpret, results. 

In §|| we briefly describe the samples we have used 
in this new analysis, including the new 7C Redshift 
Survey, which allow us to avoid problems inherent 
in the sample selection processes. In §[| we describe 
the principles of how we use these data, the types of 
objects included in the analyses, and how these prin- 
ciples were put into practice. We then examine in §^ 
the dominant trends to emerge in the various source 
properties. In §|| we present a parameterisation of our 
results for ease of comparison with previous results. 
In §^ we then review and briefly discuss results from 
the past which have hinted at the results we find. In 
order to interpret these results, we develop our model 
of radio sources in §||; within this section, we sum- 
marise in § |S.l| basic properties of radio sources which 
are widely believed to be correct and discuss what 
may be gleaned about the environments in which ra- 
dio sources reside in 58. 2. We then consider the ex- 



pansion of the source a nd o bservational constraints 
on expansion speeds in §3.3. A model for the lumi- 



nosity evolution for a given source is considered in 
§8.4; in this section we present a model for the role 
played by the hotspot in governing the luminosity of 
radio source lobes. After considering in § 3.5 the de- 
velopment of individual radio sources we describe in 
§ |8.6| how to correctly sample radio sources born at 
successive cosmic times in order to reproduce simu- 
lated versions of the data from our complete samples. 
In §8.7 we present simulated versions of our complete 
samples using the modelling described in previous sec- 
tions. In §^| we discuss what has been learned about 
the physical mechanisms by which radio sources are 
governed from our data and simulations, and briefly 
consider some implications for cosmological studies. 
We summarise the salient conclusions of our study in 

We assume that H D = 50 km s _1 Mpc -1 and, 



except where specifically stated otherwise, that Om = 
1 and fl\ = 0. We use the convention for spectral 
index (a) that S v cx v~ a , where S v is the flux density 
at frequency v. 

2. Sampling and selection 

2.1. The ineluctable biases in complete sam- 
ples 

The simplest sample-type, referred to here and in 
companion papers as a 'complete sample' has as its 
selection criteria only a single flux-limit (determined 
ultimately by the finite sensitivity of the survey tele- 
scope) in a chosen observing band and a chosen sky 
area. However, for any single flux-limited sample cho- 
sen in this way there will be an inevitable and tight 
correlation between luminosity and redshift. This 
arises^ because for a fixed flux-limit, only powerful 
sources can be detected out to great distances. It is 
therefore impossible to determine for a single flux- 
limited sample whether the dependence of a given 
source property, such as linear size, is primarily on 
redshift or on luminosity. Figure [j] shows the coverage 
of the luminosity-redshift (P~z) plane by particular 
complete samples selected at different flux limits. 

It is clear that just one such sample neither achieves 
a good range in luminosity out to high redshift nor a 
good range in redshift at low luminosity. The use of 
only one sample severely undermines any study per- 
taining to the evolution of radio sources. Extensive 
coverage of the luminosity-redshift plane is essential 
if we are to decouple the luminosity /redshift depen- 
dence of a given source property and this is only 
achievable, given restricted telescope time, using a 
number of complete samples selected at increasingly 
fainter flux- limits. 

2.2. The choice of selection frequency 

Samples selected at high frequency (> 1 GHz) come 
with complications: at these frequencies the cores of 
radio sources can be a significant contribution to the 
total flux density. Such cores are thought to be in- 
creasingly Doppler boosted as the angle between the 
jet axis and the line-of-sight is decreased (see e.g. 
Vermeulen & Cohen 1994). This means that ob- 
jects at a small viewing angle are preferentially in- 
cluded in the sample, since their total flux measure- 



1 For further discussion of the details 
tribute to this close relationship see : 
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Fig. 1. — A schematic illustration of the coverage of 
the luminosity (at 151 MHz rest-frame) versus red- 
shift (P-z) plane from a single sample selecting all 
sources brighter than 12 Jy in a certain sky area 
(coloured mid-grey), a single sample of sources with 
flux densities between 2 and 4 Jy (coloured light- 
grey) and a single sample which has selected sources 
brighter than 0.5 Jy in a much smaller area of sky 
(coloured dark-grey). Each coloured region contains 
a similar number of sources. The solid black lines 
in each case correspond to the minimum luminosity 
which a source at that redshift may have, and yet still 
be above the minimum flux-limit of the sample assum- 
ing source spectral indices are 0.5; other lincstylcs in- 
dicate the influence of other spectral indices on the 
shape of the flux-limit on the P-z plane. The up- 
per edge of the mid-grey region arises because of a 
maximum observed luminosity in radio sources. The 
upper-edge of the dark-grey region arises because in 
its small sky area the probability of finding more lu- 
minous sources is low. 



ment is promoted over the sample flux-limit because 
of their orientation, rather than their intrinsic power. 
In addition to having lower power in extended compo- 
nents, these objects will have shorter projected angu- 
lar sizes. Moreover, the GHz part of the spectrum of 
radio sources is reduced due to synchrotron losses as 
the source radiates, adiabatic expansion losses as the 
source grows and expands and, with increasing impor- 
tance with increasing redshift, due to inverse Comp- 
ton losses from the equivalent magnetic field of the 



cosmic microwave background radiation. This is fur- 
ther compounded because of the increasing difference 
between emitted frequency and observed frequency as 
redshift increases. Thus observing at a high frequency 
(> 1 GHz) will result in a vast loss in the number of 
high redshift sources which make it above the sam- 
ple limit, since inverse Compton losses increase dra- 
matically with redshift [the equivalent magnetic field 
due to the cosmic microwave background (CMB) is 
given by 0.318(1 + z) 2 nT]. Inverse Compton losses 
have been invoked to explain the X-ray flux from the 
lobes of the nearby radio galaxy Fornax A (Kaneda 
et al. 1995) and also as an explanation of the high X- 
ray luminosity of at least some distant powerful radio 
galaxies (Brunetti, Setti and Comastri 1997). Sources 
with large linear sizes at intermediate and high red- 
shifts are most unlikely to be found in high frequency 
surveys since larger sources tend to have steeper spec- 
tral indices (see §£^4). [Thus high frequency surveys 
are a poor probe of evolution in the high redshift uni- 
verse, except for the Doppler boosted population.] 

Surveys carried out at lower frequencies (a decade 
lower ~ 0.1 GHz) detect emission largely from that 
part of the spectrum which is least subject to syn- 
chrotron/inverse Compton losses and orientation bi- 
ases. They thus inform on the intrinsic energy spec- 
trum injected into the radio lobes in a way which is 
less affected by the history of that particular source. 
There is thus a strong argument for using complete 
samples detected at such frequencies. However, in 
making a radio survey at arbitrarily low frequency 
one encounters different problems: there are three 
mechanisms by which a survey conducted at ~ 10 
MHz would fail to detect sources which surveys at 
higher frequency would detect. One of these mecha- 
nisms is synchrotron self-absorption; this arises when 
the intensity of synchrotron radiation within a source 
becomes sufficiently high that re-absorption of the ra- 
diation by the synchrotron emitting particles is im- 
portant. Another mechanism is free-free absorption: 
absorption of the synchrotron radiation by thermal 
material intervening along the line-of-sight, possibly 
in the immediate vicinity of the synchrotron source 
- e.g. by its own cluster environment; an example 
of this mechanism in 3C 295 is discussed by Taylor & 
Perley (1992). 

The third mechanism arises because of a low- 
energy cut-off to the relativistic particles which give 
rise to the synchrotron emission. Arguments that 
there is a departure from a power-law in the electron- 
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energy distribution at low energies were advanced by 
Leahy, Muxlow & Stephens (1989) who made obser- 
vations of 3C radio sources at 151 MHz with MERLIN 
and at 1.5 GHz with the VLA and found spectral flat- 
tening in the hotspots of three of the powerful sources 
at low frequency. 

An extreme example of sources which are rarely 
detected in low-frequency surveys are the Giga Hertz 
Peaked Spectrum (GPS) population [O'Dea, Baum 
& Stanghellini (1991)]. It is the turnover in the spec- 
trum of these sources which gives rise to the absence 
of these sources in low-frequency surveys, and their 
presence in high-frequency surveys. 

The lower the frequency of the detection survey, 
the greater is the chance of losing small sources 
via one of these mechanisms. However, for a given 
turnover frequency, this chance decreases with in- 
creasing redshift. A sample selected at a fairly low 
frequency (such as tens of MHz) might miss small low 
redshift sources which would be detected at higher 
redshift because the observing frequency is likely to 
be less than the turnover frequency. In order to quan- 
tify the significance of this effect in samples selected at 
a frequency of 38 MHz, we considered that complete 
sa mple which is a subset of the 3C sample (described 
in § |3.l| ) defined by declination, 5 > 60°. We then 
considered which of the sources in this sample of 29 
objects would not have been included in a survey se- 
lected in the same sky area at the lower frequency of 
38 MHz, due to spectral turnover. 

The limiting flux density of the 8C survey carried 
out at 38 MHz ]Rccs (1990) ; catalogue revised by 



Hales et al. (1995)| is 1 Jy over most of the survey 



area. We set a threshold of 25 Jy for our 38-MHz 
selected sample since this is the flux density which an 
object on the 3C flux-limit at 178 MHz would have 
at 38 MHz if it had a spectral index of 0.6 — i.e. a 
spectral index which can plausibly arise without any 
recourse to one of the spectral turnover mechanisms 
discussed above. 

Three bona fide sources from a sample of 29 ob- 
jects would be lost because of spectral turnover; rele- 
vant properties of these objects are shown in Table [l]. 
Thus the deficit of small objects (none of which are 
at high redshift) missing from a sample selected at 38 
MHz rather than a somewhat higher frequency e.g. 
178 MHz is at the ~ 10% level. We therefore note 
that stronger linear size evolution might in principle 
be found from surveys at 38 MHz because of the ab- 
sence of small (low redshift) sources. 



Observing at a higher frequency will tend to detect 
a higher number of small low-redshift sources which 
have a spectral turnover, since a higher observing fre- 
quency is more likely to be above typical turnover fre- 
quencies. Even with this, we suggest that a high fre- 
quency like 408 MHz might measure spuriously high 
linear size evolution, [since a survey at 408 MHz se- 
lects objects at z ~ 2 at rest- frame 1 GHz] , since large 
and small low redshift sources will still make it into 
the sample but the only high redshift sources which 
make it into the sample are very short ones which 
are too young to have suffered catastrophic inverse 
Compton losses to their lobes. 

To investigate the effects of increasing the survey 
frequency on the selection process, we considered the 
differences between the 6C sample selected at 151 
MHz (described in §^2) and a sample selected over 
the same region of sky but at 408 MHz. We were 
able to do this using the sample of Allington-Smith 
(1982) which is roughly the same region of sky as the 
6C sample but selected at 408 MHz to be brighter 
than 1 Jy as measured in the Bologna (B2) survey. 
Those sources which are only selected at 408 MHz in 
this way are listed in Table |[ We have redshift infor- 
mation for all but one (B2 0927+35) of these objects 
hence we could overlay the points from this sample 
on our various plots of the four different source prop- 
erties (luminosity, spectral index, redshift and linear 
size; an example is presented in Figure^). We found 
that while selection of this sample at the higher fre- 
quency does result in the inclusion of more objects 
which tend to have fairly small projected linear sizes, 
these are over a range in redshift and we discern no 
redshift systematic bias implied by the properties of 
the sources included in the sample selected at 408 
MHz rather than 151 MHz, at these flux-limits. 

2.3. The use of one selection frequency 



Neeser et al. (1995) pointed out that comparing 
samples with very different selection frequencies will 
find a spurious evolution in linear size, if a faint high 
frequency (e.g . ~ 1 GHz) sample with the biases de- 
scribed in §2^ were to be compared with a bright low 
frequency (e.g. 0.1 GHz) sample without these biases. 
A study of the evolution of radio sources based on 
samples selected at different frequencies is not help- 
ful because it amounts to a comparison between the 
same types of objects at different parts of their radio 
spectra. 
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2.4. Source type 

It should be clear that samples which are incom- 
plete or inhomogeneous should not be used for a 
study pertaining to evolutionary properties because of 
the subtle and unsubtle biases which will be present. 
Even "well-defined" samples are not witho ut their 
pr oblems. For example, past studies [e.g 



Barthcl 



& J yliley (1988 )| have used samples comprised only of 
steep spectrum quasars. The motivation for such an 
approach is understandable: the spectroscopy is rela- 
tively easy hence good coverage out to high redshift 
may be achieved. However, if we believe the unifi- 
cation by orientation model [ [Bchcucr (1987)| ; [Barthcl 
(19 |39)| ] at some level (and thus that for quasars the 
jet axis is at a small angle to the line of sight), then 
the projected angular sizes which we measure will be 
heavily foreshortened. One cannot begin to depro- 
ject the sizes unless one has very large numbers of 
sources and a good idea what the mean and spread in 
the transition angle are. This problem is potentially 
compounded if quasars are used which are not from 
complete samples but have optical selection criteria 
as well, since this might preferentially bias towards 
boosted objects which are pointing towards us (see 
Willott et al. 1998a). 

In principle therefore, studies which have used only 
radio galaxies, which are presumed to lie closer to 
the plane of the sky and therefore will have a smaller 
difference between the projected size we measure and 
the deprojected size we would wish to measure, might 
be a more reliable indicator of the true dependences. 
In fact we chose to use all radio galaxies and quasars 
in our samples which fulfilled the criteria discussed in 
§^|. Such a selection criterion should ensure a random 
source orientation, allowing statistical deprojection of 
the linear size information. 

2.5. Spectroscopy 

Many previous studies of samples of radio sources 
have required the use of redshift estimators for a sig- 
nificant fraction of the sample (e.g. Dunlop & Pea- 
cock 1990). The most common estimator employs the 
fC-band (2.2 //m) magnitude. Use of this estimator 
is motivated by the small spread in the near-infrared 
Hubble Diagram (K-z relation) for radio galaxies in 
the 3C sample (Lilly & Longair 1984; Best, Longair 
and Rottgering 1998). It is now known that sim- 
ple extrapolation of the 3C relation to higher red- 
shifts and/or to fainter radio samples (which contain 



fewer radio luminous sources at any given redshift) 
leads to both large uncertainties in the estimated red- 
shifts, and to systematic biases. At redshifts of z > 2 
the K— magnitude becomes contaminated by lumi- 
nous emission lines, and other sources of non-stellar 
radiation [Eales & Rawlings (1993, 1996); Simpson, 
Rawlings & Lacy (1998)]. This is at least partially 
responsible for a scatter of 3 magnitudes in the K— 
z diagram at z > 2, and the little, if any, correla- 
tion between the variables (Eales et al. 1993). This 
large scatter may also be influenced by a genuine in- 
crease in the spread of stellar luminosities at early 
cosmic epochs. Comparison of the K-z relations of 
3C and 6C samples (Eales et al. 1997) reveals a signif- 
icant difference in the sense that 3C radio galaxies are 
about one magnitude brighter than the 6C galaxies at 
ZRil. Further studies (Best et al. 1998; Roche, Eales 
& Rawlings 1998) suggest that this effect is caused 
by a (positive) correlation between radio luminosity 
and the stellar luminosity of the host galaxy which is 
not yet well understood. To summarise, although it 
may now be possible to construct more reliable single- 
colour redshift estimators using if-band imaging of 
the combined 3C/6C/7C dataset, the use of almost 
exclusively spectroscopic redshifts has been an essen- 
tial ingredient of the work described in this paper. In 
the very few cases where we have been unable to se- 
cure spectroscopic redshifts (§3.2 and §[Q|), we have 
estimated redshifts by obtaining multi-colour (optical 
and near-infrared) information and fitting spectral en- 
ergy distributions (SEDs) with template galaxy spec- 
tra (e.g. Dunlop & Peacock 1993). This method is 
likely to be robust whenever non-stellar contributions 
to the SED can be safely neglected. Evidence that 
this is the case comes from the lack of any detectable 
emission lines in both optical and near-infrared spec- 
troscopy of the relevant objects. 

3. Samples used 
3.1. 3C 

The complete sample is that of the revised 3CR 
catalogue of Laing, Riley & Longair (1983), which 
contains 173 objects in an area of sky of 4.239 sr 
with 5 > 10° and \b\ > 10° each with a flux den- 
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sityQ greater than 10.9 Jy at 178 MHz|. We excluded 
3C231 (M82) since its radio emission is due to star- 
formation (Condon et al. 1990) rather than an ac- 
tive galactic nucleus. We also excluded another two 
objects, 3C345 and 3C 454.3, because their emission 
at 178 MHz only exceeds the flux-limit because of 
Doppler boosting of their core emission. The final 
redshift of this sample, for 4C 13.66, was obtained by 



Ra- dings et al. (1996b); the penultimate redshift, for 
4C 74.16, is 0.568 (R.A. Laing, priv. comm.). The 
spectroscopic redshift completeness for this sample is 
thus 100%. 

3.2. 6C 

The complete sample we used is based on Eales' 
(1985) sample, but re-selected from the final versions 
of parts II and VI of the 6C survey at 151 MHz [Hales 
et al. (1988) and Hales et al. (1993) respectively]. This 
revised sample is for objects with 2.0 < S151 < 3.93 
Jy (where S151 is the flux density in Jy at 151 MHz), 
with 08 20 30 < RA (B1950) < 13 01 30 and also 
with +34 01 00 < <5 (B1950) < +40 00 00. The sky 
coverage of this sample is thus 0.102 sr. Details of the 
revised sample are in Blundell et al. (1998b) and red- 
shifts are in Rawlings et al. (in prep.). Structural clas- 
sifications and angular sizes were re-measured from 
maps by Law-Green et al. (1995), Naundorf et al. 
(1992), Naundorf (1992) and Tansley (1997). The re- 
vised complete sample would consist of 59 objects but 
one (6C1036+3616) is formally excluded for a reason 
which is not statistically biassed, namely because a 
bright star lies too close to its optical identification 
to measure the redshift of the radio source. Only 
one redshift for this final sample of 58 sources is not 
spectroscopic giving 98% spectroscopic redshift com- 
pleteness for the 6C sample. 

3.3. 7C 

The 7C survey was carried out with the Cambridge 
Low Frequency Synthesis Telescope (CLFST) at 151 
MHz ]McGilchrist et al. (1990)| at a resolution of 
70 x 70 cosec(£>ec) arcsec^. in this paper we use 
the 7C-I and 7C-II parts of the 7C Redshift Sur- 
vey, which were chosen by Rossitter (1987) to over- 



2 For the analysis described in this paper, flux densities quoted 
by Laing et al. (1983) for all objects in this sample were mul 



i plied by 1.09 to put them on the same flux scale as Baars et 



0. (1977) 



lap with the 5C6 and 5C7 fields of Pearson & Kus 
(1978), respectively. The 7C 151-MHz data were sub- 
sequently re-analysed (and in the case of the 5C6 field, 
re-observed) by Julia Riley. VLA data, where avail- 
able, were retrieved from the archive (by kind per- 
mission of Barry Clark), or were recently obtained 
by us: full details of all these radio observations and 



the sample selection are described by Blundell et al 
(1998b)| . The 7C-I region covers an area of 0.0061 



sr and contains 37 sources with 6*151 > 0.51 Jy. The 
7C— II region covers an area of 0.0069 sr and contains 
40 sources with Si 51 > 0.48 Jy. All the sources in the 
final sample have now been mapped with the VLA 
giving high- resolution radio maps. We used these 
maps to deduce the structural information and an- 
gular size of each object in the 7C sample. Spectro- 
scopic redshifts have been obtained for 71 members of 
the 7C complete sample, thus spectroscopic redshift 
completeness is 92% (Willott et al. in prep.). For the 
remaining objects, their redshifts are estimated from 
4 or 5 optical/near- infra red magnitudes as described 
in Rawlings et al. (in prep.). Only one source in the 
small sky area of the 7C sample (3C200) is common 
to the 3C sample, and this source is just included 
as a member of 3C for the purposes of the analysis 
presented in this paper. 

We excluded one object, 5C7.230 (also known as 
7CB082134.3+244829), from the sample because its 
emission at 151 MHz only exceeds the flux-limit be- 
cause of Doppler boosting of its core emission (Blun- 
dell et al. 1998b). 

4. Classifications & measurement techniques 
4.1. Source structure. . . 

Neeser et al. (1995) demonstrated the perils of try- 
ing to analyse the changes in maximum angular extent 



of the 'edge-darkened' FRI-class [Fanaroff & Riley 



(1974)| radio sources with increasing redshift and at- 



tributed the inclusion of such objects to the extremely 
strong linear size evolution with redshift which has 



been reported in the past (e.g. Port, Katgert fc Wind 



horst (1987) ). To avoid spurious correlations and to 



3 This roughly corresponds to a flux-limit of 12 Jy at 151 MHz. 



focus our analysis on the evolution of one physically 
distinct class of object out to very high redshift, i.e., 
to study the same physical phenomenon at different 
redshifts in order to compare like with like, we inves- 
tigated the radio spectra and structure of the objects 
in our samples. One analysis which is part of this 
paper — the linear size evolution of classical double 
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radio sources — amounts to considering how the sep- 
aration of working surfaces at the ends of oppositely 
directed radio jets depends on luminosity and/or red- 
shift. For this reason we excluded from our analysis 
radio sources of the FRI structural type. We regard 
objects of this morphological type as the inevitable 
manifestation of low jet-power objects, where the jet- 
power is insufficient to form a supersonic shock, given 
its environment, where the jet impinges on the inter- 
galactic medium [see e.g. Falle (1991)]. The 6C sam- 
ple and 7C sample selected to be a factor of 6 and a 
factor of nearly 25 respectively fainter than 3C, in fact 
contain relatively few objects of FRI structure. This 
arises because in the relatively small sky areas of 6C 
and 7C few nearby bright but low-luminosity objects 
are found. Note that luminosity is not a fundamental 
quantity in these objects in the way that the power 
transported by the jet is. Distant low-luminosity ob- 
jects which are found will not be the exact counter- 
parts of low-redshift objects of the same luminosity: 
they must have higher jet power to compensate for 
the fact that for an object of redshift z, inverse Comp- 
ton losses increase as (1 + z) A . Put another way, an 
object with a given jet power observed in an environ- 
ment with a given density profile at high redshift will 
have a much lower luminosity at, say, 80% of the way 
through its life than a lower redshift object with oth- 
erwise the same parameters at the same point. [Note 
that this is not an observational error: the selection 
is made at a sufficiently low frequency and resolution 
that there is no question of missing sources because 
of undersampling smooth extended flux.] 

For this reason we do not base the criteria of which 
objects should be included in our analysis on their 
luminosity, but on their structure. 

How then to identify whether a given source is type 
FRI or FRII? We deemed that an object was a clas- 
sical double (FRII) if at least one hotspot was at the 
leading edge of a working surface. In practice this 
meant a bright, compact feature with a sharp discon- 
tinuity in surface brightness was at the leading edge 
of emission expanding outwards from the core or loca- 
tion of the optical identification. We did not require 
that this was the case on both sides of the core [see e.g., 
5C6.62 also known as 7C B021029. 7+325402 (Blun- 
dcll et al. 1998b)] to allow for the fact that although a 
working surface can exist, it can have a slightly weird 
appearance due to e.g., Doppler suppression, but the 
presence of one such hotspot is a good indicator that 
in the environment in which a given source is located, 



the jet-power is sufficiently high that a classical dou- 
ble radio source will be manifested. 

In some objects there was clear evidence of emis- 
sion beyond a possible working surface, which had 
decreasing surface brightness with increasing distance 
from the core. We considered these to be a somewhat 
different physical phenomenon from the 'working sur- 
face + backflow towards the core' type object, and 
called them FRIs and excluded them from our anal- 
yses. For those intermediate cases where we could 
not be totally sure that an object should be assigned 
to the FRII category, we considered by examination 
of their integrated radio spectra (see §4.2) whether 
there was any evidence of trailing outflow with in- 
creasing distance from the core, as might be evidenced 
by a greater flux from a survey measurement (e.g., 
from VLA D-array or single dish measurements) than 
from our higher resolution (VLA A- or B-array) mea- 
surements. Thus, objects which do not apparently 
contain any clearly defined working surface, but in 
which we believe all the flux is contained within a 
well-defined measureable extent, are classed as having 
defined angular sizes (DAs). We believe that most, if 
not all, the angular sizes measured for these objects 
do correspond to the separation of working surfaces 
at the ends of oppositely directed radio jets, though 
we cannot eliminate the possibility that some might 
be FRIs. 

Our angular size measurements were our best at- 
tempt at summing the separations of the two working 
surfaces from the core; in the absence of a radio core 
we took our best estimate of the separations of the 
working surfaces from one another. 



4.2. 



, and Spectra 



Fitting of the radio spectra for the members of all 
three samples was made possible by consideration of 
data from many radio surveys^, and in the case of the 
7C sample in combination with data from our own 



VLA maps as tabulated in Blundcll et al. (1998b) 



We fitted the flux density data using a Bayesian 
polynomial regression analysis which assessed the 
posterior probability density function (pdf) for the 
required order of a polynomial fit (e.g. Sivia 1996); we 



4 In addition to the original survey data at 151 or 178 MHz, 
we obtained from the literature flux densities at 365 MHz, 408 
MHz, 1.4 GHz and 4.86 GHz, and in the case of 3C we also 
obtained flux densities at 10, 22, 38 and 85 MHz, chiefly from 
Laing & Peacock (1980). 
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chose x = log 10 (^/MHz) as the independent variable 
and expanded the dependent variable y = log 10 &v = 
12^=o a rZ r - I n a number of cases the peak in the 
PDF led us to prefer a first-order polynomial, and 
the gradient of the fit (<ii) was multiplied by —1 to 
give the spectral index. In many cases the PDF led 
us to prefer a second-order fit, implying significant 
curvature in the radio spectrum: in these cases we 
calculated the curvature (3 = — 2a2, and hence the 
spectral index at (rest frame) 151 MHz which equals 
— a\ + /?log 10 (151/(1. + z)). In no case was a higher 
order fit preferred. 

With a fitted functional form for the spectrum 
of each object, we interpolated (in the case of 3C) 
or extrapolated (in the case of 6C and 7C) to de- 
rive the flux densities at rest-frame 151 MHz and to- 
gether with the redshifts we derived the luminosities 
in WHz -1 sr^ 1 in a number of different assumed cos- 
mologies, according to the formulae presented in Car- 
roll, Press & Turner (1992). 

5. Results 



We derived all plots in the following three assumed 
cosmologies: (1) 12m = and Qa = 0, (2) £1m = 1 and 
fl\ — and (3) = 0.1 and Qa = 0.9; similar de- 
pendences were seen in each case. We also considered 
separately just the class of certain FRIIs and also the 
combined set of FRIIs and DAs (see §4.1). No signif- 
icant difference was found. The plots we present here 
are for the combined set of FR lis and DAs and, ex- 
cept for Figure]^, an assumed cosmology with = 1 
and £!a = but the statistical analysis presents the 
numbers for all three cosmologies we considered (see 



5.1. Plots of source dependences 

Figure ^| shows the luminosities vs projected linear 
sizes (the P-D plane) of the radio sources. No strik- 
ing dependence of the luminosity of an object on its 
linear size is seen, though there is a deficit of small, 
low-power objects in the lower-left quadrant of the 
P-D plane, and also of large, high-power objects in 
the upper-right quadrant of the P-D plane. It can 
be seen that the region most densely occupied by the 
solid circles (corresponding to objects with z > 2) is 
displaced to the left [smaller projected linear sizes] of 
the region occupied by the main density of plus signs 
(corresponding to objects with z < 0.8), indicating a 
distinct anti-correlation of linear-size with redshift. 




0.1 1 10 100 1000 

Linear size (kpe) 



10000 



Fig. 2. — This graph shows for all members of the 
complete samples included in our analysis their lu- 
minosities plotted against their projected linear sizes. 
These properties are calculated assuming S1m = 1 and 
f?A = 0. Four different redshift bins are indicated by 
symbol type. 

Figure || shows more directly that sources at higher 
redshift tend to have shorter projected linear sizes. 
This trend is seen throughout the redshift range, and 
is seen in each of our three considered cosmologies. 
The symbol type in this figure delimits the luminos- 
ity bin an object is in. In the middle plot where the 
assumed cosmology is Qm = 1 and CI a — 0, the most 
powerful sources are not as powerful as in the cos- 
mologies with low CIm- This plot demonstrates that 
in breaking the degeneracy between luminosity and 
redshift, above a redshift of 2 the tendency for the 
sources to be smaller seems to be independent of the 
fact that the majority of these sources are not in the 
highest power bin. 

The operational definition of spectral index used in 
the past has been to approximate the radio spectrum 
of a source by a single power law and to derive a as 
— dlogSu/dlogv, invariably in the observed frame, 
using flux density measurements at two frequencies. 
When this usage of spectral index is employed, there 
is a significant spectral index dependence on redshift 
(see Figure ||) which has been exploited (Chambers, 
Miley & van Breugel 1990; Blundell ct al 1998a) to 
find the most distant radio galaxies. This precise 
technique was used to find the most distant radio 
galaxy currently known (Rawlings et al. 1996a). Once 
the curvature and redshifting of the spectra ("K- 
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Fig. 3. — These graphs show for all members of the 
complete samples included in our analysis their pro- 
jected linear sizes against their redshifts, with these 
properties calculated for three different assumed cos- 
mologies. The luminosity bin an object falls into is 
denoted by symbol, with the same threshold values in 
each case. 

correction" ) are taken account of, this dependence of 
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Spectral index (5 GHz rest frame) 

Fig. 4. — These graphs show for all members of 
the complete samples included in our analysis red- 
shift against spectral index measured in two different 
ways. The top panel shows the spectral indices before 
fitting for spectral curvature and i^-correction and 
assuming a single power-law between 151 MHz and 
1.4 GHz while the lower panel shows spectral indices 
evaluated at 5 GHz in the rest frame of each source. 
Source luminosities are indicated by symbol-type. 



spectral index on redshift is dramatically different, in 
a way which depends on the frequency at which the 
spectral index is evaluated. When the spectral in- 
dex is evaluated at 5 GHz rest-frame, then a rather 
strong dependence on redshift is seen (see Figure 
This arises because the GHz part of the spectrum is 
where inverse Compton losses (as well as other loss 
mechanisms) are most manifest. If spectral indices 
are evaluated in the rest frame at a frequency of a 
few hundred MHz, this strong dependence on redshift 
diminishes (see Figure ||) . 

Figure || shows a plot of the luminosity of each 
object against its spectral index, each evaluated at 
rest-frame 151 MHz. There is a clear trend for the 
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Fig. 5. — These graphs show for all members of the 
complete samples included in our analysis their red- 
shifts against spectral indices measured in two differ- 
ent ways, in both cases fitting for spectral curvature 
and X-correction. The upper panel has the spectral 
indices evaluated at 151 MHz in the rest frame while 
the lower panel shows the spectral indices evaluated 
at 300 MHz in the rest frame of each source. Source 
luminosities are indicated by symbol-type. 

most powerful objects to have the highest spectral 
indices. 

At a frequency as low as rest-frame 151 MHz, the 
spectral index is informing us on the energy index of 
the synchrotron particles as injected into the lobes. 
We hereafter refer to this power-law exponent of the 
energy distribution of particles injected into the lobes 
as the 'injection index', reserving the term 'accelera- 
tion index ' for the exponent of the energy distribution 
to which particles in the hotspot may accelerated [and 
described by a model such as diffusive shock acceler- 
ation e.g. Bell (1987)]. We explore in §[0 how this 



29 - 

28 - 

27 - 

26 - 

25 - 
24 _ 




• z a 2 

o 2 > z 1 1.5 

X 1.5 > z * 0.8 

+ 0.8 > z 



injection index may depend in a simple way on the 
intrinsic jet-power of an object and hence how such 



-0.5 0.0 0.5 1.0 1.5 

Spectral index (151 MHz rest frame) 

Fig. 6. — This graph shows for all members of the 
complete samples included in our analysis their lu- 
minosities plotted against their spectral indices each 
evaluated at rest-frame 151 MHz. The redshift bins 
are indicated by symbol type. 

an injection-index-luminosity dependence can arise. 

A neat manifestation of the P-a correlation is seen 
in the luminosity-redshift plane. Figure shows lines 
which trace the flux-limits for 3C, 6C and 7C, i.e. 
12, 2 and 0.5 Jy respectively. The solid black lines 
in each case correspond to the minimum luminosity 
which a source at a given redshift may have, and yet 
still be included in a given sample if its spectral index 
is 0.5; other line-styles indicate the influence of other 
spectral indices on the shape of the flux-limit on the 
P—z plane. It can be seen that those 3C sources with 
luminosities above ~ 10 27 WHz -1 sr _1 all avoid the 
region between the flux-limit line for a spectral index 
of 0.5 and that for a spectral index of 0.8. For those 
3C sources which have rather lower luminosities, it is 
impossible to distinguish any preferred flux-limit line. 
In this lower luminosity regime the fainter samples 
do not show any such avoidance of the flux-limit line 
given by a spectral index of 0.5. 

Inspection of the D—a plane depicted in Figure || 
shows a correlation of projected linear size with spec- 
tral index, although it also suggests that sources with 
even quite large linear sizes may not necessarily have 
very steep spectra at rest-frame 151 MHz. It is pos- 
sible that some of these spectra are flattened by free- 
free absorption from intervening material. Taylor & 
Perley (1992) present a study of 3C295 which has a 
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Fig. 7. — These graphs show for all members of the 
complete samples included in our analysis their lumi- 
nosities at rest-frame 151 MHz against redshifts. The 
top panel shows the P—z plane for the 3C and 6C 
samples; the lower panel shows the P—z plane for the 
7C sample. 



(complex) turnover in the spectrum below 100 MHz; 
they attribute this to thermal absorption from multi- 
component ISM gas at the cluster. 

5.2. Disentangling the different dependences: 
statistics 

We display in Tables || and || the results of the 
statistical analysis using the 4-way Spearman partial- 
rank correlation coefficient described by Macklin (1982) 
to indicate the true dependences of the different 
source properties. In these Tables, r a ^ cc i 
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Fig. 8. — This is a plot of projected linear size against 
spectral index evaluated at 151 MHz rest-frame for 
all members of the complete samples. We also over- 
lay those sources which would have been members of 
the 6C complete sample had the selection frequency 
been 408 MHz [taken from Allington-Smith (1982); 
see §|2.2[]. The redshifts of the other sources are indi- 
cated by symbol type. 

lation coefficient between o and b when c and d are 
held constant. We repeat this analysis for the same 
three assumed cosmologies, (i) f^M = and = 0, 
(it) Qm = 1 and f^A = and (Hi) Qjyr — 0.1 and 
f^A = 0.9. The dependences largely persist in all the 
cosmologies; we quantify and discuss the meaning of 
the dependence on cosmology of the strength of the 
anti-correlation of linear size with redshift in a forth- 
coming paper. 



5.3. A cautionary tale: 
limited sample 



using just one flux- 



is the corre- 



To demonstrate how the tight correlation between 
luminosity and redshift in a single flux-limited sam- 
ple can obfuscate the true dependences between these 
source properties and linear size and spectral index, 
we show in Table || the results of the statistical analy- 
sis carried out above on only the 3C sample (§|3.l|). Al- 
though the relationship between linear size and spec- 
tral index is still discernible, (indicating that it is 
a real effect independent of luminosity or redshift) 
the extremely strong correlation between luminosity 
and redshift swamps the other dependences which do 
emerge when a number of complete samples with dif- 
ferent flux-limits are taken together. 
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6. Parameterisation of the linear sizes 

Although we later demonstrate that there is no di- 
rect single physical mechanism by which linear size 
depends on redshift, for comparison of our data with 
work previously done by others we calculate the para- 
metric dependence of the linear size on redshift and 
on luminosity with the form: 



D oc (1 + z)~ n P r ' 



(1) 



Following Neeser et al. (1995), values of the ex- 
ponent n are found by multiplying individual source 



siztis by (i+z) trntn Qie euiielauuii eueilicieiil (ueiice 
its associated significance) gucs Lu zciu. — A funnal 

error-bar is derived by considering the values of n 
at which the significance of the correlation coefficient 
equals unity. The exponent m was found in a simi- 
lar manner. In Table | we list the various exponents 
thus derived. The notation we use is that n[r a f,| c( i] (or 
m [ r ab\cd)) is the exponent required to obtain a corre- 
lation coefficient of zero between quantities a and b 
when c and d are held constant. For comparison with 
previous studies, which have not included the spec- 
tral index dependence, we also present 7"i[r o;) | c ] and 
m[r a h| c ] which are the exponents required in equa- 
tion p] to make the 3-way correlation coefficient be- 
tween a and b equal zero, when c is held constant. 



7. Comparison with previous work 
7.1. The variation of linear size 

It has never been in dispute that the linear sizes of 
classical double radio sources are larger in the local 
Universe, although this has often come from studies 
restricted to very low redshifts. Indeed, this together 
with the fact that the angular sizes of classical double 
radio sources are relatively easy to measure led some 
to hope that it might be possible to use these objects 
as ' standard rods' to co nstrain cosmic geometry (see 
e.g. Ubachukwu (1995) ), if it could be demonstrated 
that the linear size of an object is not influenced by 
its redshift. However, past reports of the true depen- 
dence of the linear sizes of radio sources have been 
equivocal as a brief survey of the literature indicates: 



Kapahi (1987) found a strong evolution in linear 
sizes of the form D cx (l + z) _3±0 5 for a cosmol- 
ogy with = 1 and f^A = 0, but at constant 
redshift he found that D oc p°- 3 ± 01 . This was 



for a set of FR II radio galaxies from four sam- 
ples: (i) The BDFL sample (Bridle et al. 1972) 
selected at 1.4 GHz; these galaxies had almost 
complete spectroscopic information and have 
0.08 < z < 0.2. (ii) The MC-I and MC-II 
samples (Machalski & Condon 1983; Machal- 
ski & Condon 1985) selected at 1.4 GHz. Out 
of these two samples only 4 galaxies had spec- 
troscopic redshifts; redshift estimates from R- 
band magnitudes were used and redshift limits 
inferred for the two samples were 0.4 and 0.6 
respectively. (Hi) The LBDS sample which was 



selected at 1.4 GHz by Windhorst, van Heerde 



Katgert (1 



For the analysis Kapahi pre- 
sented, all sources in this sample were deemed 
to have a redshift of 1.25. (iv) Galaxies from 
the revised 3CR sample which was selected at 
178 MHz by Laing et al. (1983). 



Barthel & Miley (1988) found a weaker, and 
more explicitly uncertain, dependence of linear 
size on redshift and essentially no dependence 
on luminosity (D oc (1 + z )-i-8±i-*p-0-03±o.sj 
from a sample of 134 steep spectrum radio 



quasars from Hewitt & Burbidge (1980) known 
to exhibit radio emission. In the context of 
finding for their sample that distant quasars 
have a more bent and distorted appearance than 
nearby ones, they postulated that the ambient 
medium, interacting with the radio jets, played 
an important role in governing linear size evo- 
lution. 



• Singal (1993) found for radio galaxies in his 
chosen sample that there was a strong anti- 
correlation of linear size with redshift but a cor- 



relation with luminosity [D oc (1- 



-3.0 



p+0.35j 



while for quasars an anti-correlation of linear 
size and luminosity was found as was marginal 
evidence for a dependence of linear size on red- 
shift. His sample consisted of 789 radio sources, 
taken from 3CRR selected at 178 MHz by Laing 
et al. (1983), the 1-Jy sample selected at 408 
MHz (Allington-Smith et al. 1988; Lilly 1989), 
and the Molonglo (McCarthy et al. 1991) and 
B3 samples (Vigotti et al. 1989; McCarthy 
1991) which were selected at 408 MHz. 

• Oort, Katgert & Windhorst (1987) found that 
the dependence of linear size on luminosity and 
on redshift, assuming that Qm — 0, was de- 
scribed by: D OC (l + )-3.3±O.5pO.3±O.O5_ Then . 
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analysis was based on 163 radio galaxies from 
a number of different samples which included 
FRI type radio galaxies. 



• Nilsson et al. (1993) found from a sample of 540 
FRII double sources taken from the literature 
that for (high power) radio galaxies and quasars 
linear size and power are anti-correlated; they 
did not find any evidence for linear size evolu- 
tion. 

• The first study to use only complete samples, 
and using almost entirely spectroscop ic red- 
shifts, was that by Neeser et al. (1995)| . Using 
the rev ised 3CR an d 6C samples [Laing et al. 
(1983); |Eales (1985| , described in §|1] and §|J 
respectively, they found D cx (1 + z)~ 1,2±0 ^ for 
n M = and tt A = and D cx (1 + 2 )-i-7±o.5 
for r2jvi = 1 and Cl\ — 0. They found no in- 
trinsic correlation between size and radio lu- 
minosity. Neeser et al. (1995) suggested the 



presence of dense line-emitting gas as being re- 
sponsible for the observed linear size evolution. 
They also pointed out that such an effect could 
be achieved if radio jets were switched on for 
shorter periods of time at high redshift than at 
low redshift. 

A recent study (Buchalter et al. 1998) was based 
on the 1.4 GHz FIRST survey; this study fo- 
cussed on a sample of 103 quasars from the 



quasar catalogue of Hewitt & Burbidgc (1980) 



which had radio counterparts in the FIRST sur- 
vey. In complete contrast with Neeser et al 



(1995), they found no evidence for a dependence 
of linear size on redshift, other than that which 
arose from an anti-correlation of linear size with 
luminosity. 

It is clear that sample inhomogeneity and redshift 
estimates rather than measurements, as examined in 



detail by Neeser ct al. (1995) , have played a significant 
role in giving rise to this plethora of contradictory 
results. 

7.2. Spectral index dependences 

Almost 40 years ago, evidence for a relationship 
between spectral index and luminosity was suggested 
by Heeschen (I960) . Since that time other studies 



have confirmed his initial findings, e.g. Blumcnthal 



they were unable to identify whether the dominant 
dependence was on redshift or on radio power. 

Laing & Peacock (1980) reported that for sources 
with hotspots, spectral index and the 'degree of spec- 
tral curvature' were correlated with luminosity. Their 
study was based on two complete samples of radio 
sources: one was 3CR (Bennett 1962) selected at 178 
MHz and the other was selected at 2.7 GHz [later 
published by Peacock and Wall (1985)]. They were 
similarly unable to say whether the dependence of 
spectral index was predominantly with luminosity or 
with redshift. Krolik & Chen (1991) discussed a num- 



ber of causes of these results: (i) since a typical radio 
galaxy spectrum steepens to higher frequency, in part 
a simple redshifting effect will give rise to the observed 
correlation between redshift and spectral index. They 
also pointed out (ii) the importance of inverse Comp- 
ton losses in higher redshift sources which will give 
rise to steeper spectrum sources and (Hi) the fact 
that an increased magnetic field which will increase 
the synchrotron emission (and hence increase the like- 
lihood of a source being above a survey flux-limit) 
will also increase the rate at which spectral steepen- 
ing occurs. No if-corrections were made to Laing 
& Peacock's (1980) spectra so the correlation they 
found probably arose from some combination of all 
three of the causes discussed by Krolik & Chen (1991). 
Employing (a large proportion of) redshift estimates 
for a faint (8C-selected) sample, Lacy et al. (1993) 
showed that high-frequency (2 GHz) spectral index 
correlated more closely with redshift than with lumi- 
nosity pointing towards the dominant importance of 
inverse Compton losses on high-frequency spectra. 

8. Modelling the complete samples 

In order to gain insight into the distribution of ob- 
jects in the six planes [P-D, P-z, P-a, D-z, D-a 
and z-a] of the [PDza] parameter space occupied 
by the radio sources which are members of our com- 
plete samples, we performed Monte-Carlo simulations 
which we describe later in this section. We begin 
in §8.1 with a precis of assumptions about the basic 



properties of radio sources which bear the test of ob- 
servations and numerical simulations with particular 



reference in §3.1.3 to the distinct roles of the hotspot 
and the head in a radio source. In §8.3.2 we discuss 



the constraints we have on the way in which radio 
sources expand with time giving us the formalism for 



k Miley (1979) studied 3C and 4C sources, though 



their linear sizes. We describe in §3.4 the change of 
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the luminosity (and the radio spectra) with time, of 
classical double radio sources and produce tracks for 
individual objects in §|8.5[ We then explain in §8.£ our 



formalism for obtaining the redshift distribution, and 
then how the model properly takes into account the 
sample selection process. We discussed in §|2| how it 
was imperative that a combination of complete sam- 
ples should be used to infer the true dependences of 
radio source properties discussed in §g. A necessary 
condition for a satisfactory model of radio sources to 
explain the dependences from the combined samples 
is that it can reproduce the [PDza] parameter space 
for each single flux-limited sample individually. We 
therefore conclude this section with a presentation of 
the six planes from the [PDza] parameter space for 
each of a simulated and real dataset for the 3C and 
7C complete samples. 

8.1. Assumptions about radio sources 

8.1.1. Basic formalism for synchrotron radiation 

The synchrotron luminosity of an emitting region 
with volume V, threaded by magnetic field B is given, 
in units of WHz -1 sr -1 , after averaging over pitch 
angles by: 



1 B 2 7 3 , . 

Pu = n{"f)V, 

D7T Z[l V 



(2) 



where 7 is the Lorentz factor of the electrons radiat- 
ing at frequency v, 71(7) is defined in terms of 71.(7)^7 
as the number of electrons per unit volume per unit 
Lorentz factor between 7 and 7 + e?7, a c is the Thom- 
son scattering cross-section, /i is the permeability of 
free-space and c is the speed of light. Adopting a 
delta-function approximation for the spectral emis- 
sion from a given relativistic electron, a given fre- 
quency of emission v is related to the Lorentz factor 
of the radiating electron by: 



7: 



eB 



-2-kv 



(3) 



where m is the rest mass of an electron and e is 
the charge on an electron. Although it is the case 
that a delta function profile is not a good representa- 
tion of the power spectrum radiated by a synchrotron- 
emitting particle (particularly in the high -7 regime), 
as long as the underlying energy distribution of the 
radiating particles is smooth, then the energy emit- 
ted by the population at a particular frequency can be 



reasonably calculated using this approximation. The 
approximation will be particularly poor close to en- 
ergy cutoffs. 

8.1.2. Jet behaviour 

We assume throughout a model in which a jet em- 
anating from the central engine of an Active Galac- 
tic Nucleus (AGN), which transports a total power 
Q , remains relativistic all the way from the central 
engine to the jet-shock (see Wardle & Aaron 1997) 
and that the mechanism for confining the jet is via 
pressure balance with the lobe. We further assume 
that the jet has a density much lower than the exter- 
nal density [supporting evidence for this assumption 
comes from the simulations of Lind et al. (1989) and 
Clarke, Norman, & Burns (1989)]. In the case of a 
low density jet, simple balancing of the momentum 
flux enabled Begelman, Blandford & Rees (1984) to 
argue that the jet-shock advances at much less than 
the speed of light and consequently a large fraction of 
the bulk kinetic energy must be converted to thermal 
energy in the post-jet shock. We assume as indicated 
by the simulations of Lind et al. (1989) that the jet 
thrust is applied (even instantaneously) over an area 
significantly larger than the cross-sectional area of the 
jet itself and we identify this region with the hotspot, 
as discussed in the next section. 

8.1.3. Heads and hotspots 

We now carefully describe our usage of the term 
hotspot: we do not use the term hotspot merely to 
mean the bright emission in the outermost parts (the 
'heads') of the lobes. We wish to refine the definition 
of hotspot from that which arises purely empirically 
(see e.g., Bridle et al. 1994). We use the term to 
refer to that region in which a distinct physical phe- 
nomenon occurs, namely emission associated with the 
post-jet shock structure (i.e. the extremely high mag- 
netic field region within and just beyond the shock 
structure where a very large fraction of the bulk ki- 
netic energy from the jet is thermalised and where 
particle acceleration may occur) . A cartoon illustrat- 
ing this is shown in Figure |^. 

Although many high-resolution high-fidelity im- 
ages of hotspots reveal complex structure (Leahy et 
al. 1997, Hardcastle et al. 1997 and Black et al. 1992), 
we find no evidence in a comparison of the high- 
luminosity quasars studied by Bridle et al (1994) with 
the lower-luminosity objects studied by Leahy et al. 
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Fig. 9. — A schematic illustration of the basic fea- 
tures of a radio source which we wish to model: one 
half of a radio source is shown. Close to the outer- 
most edge of the radio source is the hotspot where the 
bulk kinetic energy of the jet is thermalised and where 
particle acceleration may take place; the hotspot feeds 
the 'head' of the source, and it is the head pressure 
which governs the rate at which the source grows into 
its environment. The decreasing depth of the grey 
colour in the lobe represents the decrease in pres- 
sure between the head of the source back towards the 
active nucleus — the pressure gradient which drives 
backflow of plasma from the head. Note that depend- 
ing on the variation of ambient density with radius, 
that the pressure may increase once more as the host 
galaxy of the AGN is approached. 

(1997) of any very strong dependence of hotspot size 
on lobe pressure; the hotspots are of order a few kpc 
in diameter in all known FRII sources. 

We also note that whatever the complexities of 
the physics by which hotspots are governed, hotspots 
do not appear to grow self-similarly with the radio 
source: i.e. the hotspot sizes do not scale strongly 
with source length. We note that what we perceive 
hotspot sizes to be will depend on the extent to which 
the emission from each hotspot is Doppler boosted. 

8.1.4- Hotspot observables 

The properties of hotspots which we are able to 
measure from high resolution maps made with MER- 
LIN or the VLA include the flux densities and the 
break frequencies in hotspots at GHz frequencies, an- 
gular sizes subtended by hotspots (which in reality 
are probably lower limits because Doppler boosting 
of the plasma emerging through the hotspots will 
tend to make hotspots appear smaller). From these 
properties we are able to estimate the magnetic field 
strengths and pressures. (Estimates of hotspot mag- 
netic fields may be made using the assumption of 



equipartition in the hotspots. X-ray detections of 
hotspots now give considerable confidence in the va- 
lidity of the assumption of equipartition magnetic 
field strengths in hotspots (values which are in the 
region of tens of nT) [see Harris, Leighly & Leahy 
1998].) 

We describe in § 8.1.5| our justification for neglect- 
ing the direct contribution of the hotspots' luminos- 
ity to the overall luminosity measured from a radio 
source at a survey frequency like 151 MHz. Nonethe- 
less the hotspot and its magnetic field play a vital role 
in governing the energy spectra of particles injected 
into the lobe, and hence the spectrum of the total 
radio luminosity we observe. 

8.1.5. The contribution of hotspots to the total lumi- 
nosity 

We contend that at low frequencies like 151 MHz, 
the luminosity from a hotspot never dominates over 
the luminosity from the head and the lobe. Naive 
application of equation ^, using typical magnetic 
fields inferred for hotspots (e.g. as derived from equa- 
tion [ll]), for an energy distribution which extends 
down to particles with Lorentz factors of ~ 1, for 
typical hotspot volumes (~ a few kpc in diameter) 
predicts for the highest jet-power objects hotspot lu- 
minosities which are higher than the observed total 
luminosities of radio sources. However, for frequen- 
cies of a few GHz, the predicted luminosities appear 
to be much more closely matched with the lumi- 
nosities measured for resolved hotspots in observa- 
tions. Leahy et al. (1989) found for a sample of 3C 
sources that the highest power sources had hotspots 
whose spectra became flatter between 151 MHz and 
1.5 GHz. (For the calculated magnetic field strengths 
the synchrotron self-absorption turnover frequencies 
for emitting volumes with sizes of a few kpc are in 
the region of tens of MHz.) They interpreted this as 
the manifestation of a low-energy cutoff in the dis- 
tribution of the radiating particles giving a deficit of 
hotspot emission at 151 MHz. In the context of an 
electron-proton jet where particle acceleration occurs 
at the hotspot such a cutoff is likely to occur when the 
Lorentz factor is of the order of the ratio of the proton 
mass to the electron mass, since for a particle to be 
efficiently accelerated, its gyro-radius must be large 
compared with the shock thickness (which is of the 
order of the gyro-radius of the thermal material, i.e. 
the protons, present). Although the precise nature of 
the low-energy cutoff depends on details such as the 
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incoming jet velocity, it is likely that the value of the 
minimum Lorentz factor does not change greatly from 
hotspot to hotspot. However, the cutoff in frequency 
at which this is manifested does vary from hotspot 
to hotspot since for a given minimum Lorentz factor 
the minimum emitted frequency is proportional to the 
magnetic field in the hotspot (equation ||) . For those 
sources whose hotspots have higher magnetic fields 
(which we will argue via equation [II] are those with 
higher jet-powers), the minimum frequency at which 
significant radiation will be observed will be higher 
than for sources with lower hotspot magnetic fields 
(and lower jet-powers). Such reasoning is consistent 
with the low contribution (0.04%) the hotspots make 
to the total 327-MHz luminosity in the archetypically 
high-Q object, Cygnus A (Carilli et al. 1991). 

8.1.6. Inverse Compton losses arising from the AGN's 
own radiation field 

We now consider whether, for simulating low radio 
frequency surveys of the type described earlier in this 
paper, it is acceptable to neglect energy losses due to 
inverse Compton scattering off the AGN's own radi- 
ation field. This mechanism was invoked by Blundcll 
& Lacy (1995) to explain the steep spectrum in the 
GHz regime of the extended radio emission around 
the radio-quiet quasar E 1821+643, and by Brunctti 
et al. (1997) to explain the soft X-ray emission de- 
tected in six 3C FRII radio galaxies. An 'equivalent 
magnetic field' may be derived from the energy den- 
sity in the radiation field; this equivalent magnetic 
field falls off as l/r where r is the distance out from 
the AGN. Significant losses arising from this mecha- 
nism will therefore be confined to some specific spatial 
extent. If this mechanism causes spectral steepening 
to be manifested at frequencies which may be red- 
shifted to be close to the survey frequency, then its 
contribution cannot be neglected; we now investigate 
whether this could be the case. 

The energy density of the AGN radiation field 
[/ rad is given by: J7 rad = Q P hot/(47rcr 2 ) where Q p hot 
is the photo-ionising luminosity. We assume that 
Qphot ~ 4>Qo where Q is the total kinetic energy 
transported by one jet and is a factor which Willott 
et al. (1998b) suggest might be as high as 20. Equat- 
ing the f7 ra d with /2fi we find, 



B eq (r) 



1 



-Mo 



27T7- 2 



(4) 



over which this magnetic field will impact on the flux 
which may be detectable by a survey, at a given dis- 
tance from the AGN. We plot in Figure [II] a thick 
solid line indicating the trajectory of a point moving 
out from the centre of the AGN at which the break 
frequency would lower to 1 GHz (this will not di- 
rectly impact on a survey at a frequency as low as 
151 MHz, until the source redshifts are greater than 
5) if the plasma bathed in this radiation field had 
resided there since the radiation field was switched 
on. [This therefore over-estimates the rate at which 
the break in the spectrum will advance out from the 
centre of the AGN.] This trajectory is calculated us- 
ing equations |l2| and ^ at each position r. We also 
plot a thin solid line which is the trajectory of one 
head advancing away from the AGN (as described by 
halving the total linear size expressed in equation ^ 
while the dashed line is the trajectory of the mid- 
point of one arm of the radio source. It is seen that 
the spatial point where the break frequency becomes 
as low as 1 GHz is only reached after the midpoint 
of the lobe has advanced beyond this. Comparison 
of the Inter-Planetary Scintillation measurements of 
the low-frequency (81.5 MHz) measurements of 3C 
objects by Purvis et al. (1987) with the flux densi- 
ties of these objects at 86 MHz of Laing & Peacock 
(1980) suggest that roughly half of the low- frequency 
emission comes from the head of the source. For the 
other half of the emission coming from the lobe, a de- 
creasing proportion comes from plasma with increas- 
ing proximity to the core. Thus, because of the speci- 
ficity of the spatial extent of the equivalent magnetic 
field arising from the quasar radiation field, those 
regions which suffer the most from this loss mecha- 
nism are those which barely make a significant con- 
tribution to the overall low-frequency flux. The im- 
portance of this effect cannot be so simply analysed 
for sources which are smaller in extent than e.g. 10 
kpc. However, such sources are known to expand 
very quickly [Owsianik & Conway (1998)] and the ra- 
diating plasma in them is continually replenished by 
freshly injected synchrotron emitting particles. Thus 
we conclude that the effects of inverse Compton losses 
off the AGN's own radiation field should not be sig- 
nificant in the selection of sources detected in surveys 
conducted in these frequency regimes, although this 
mechanism may be relevant in the question of the 
'docked tails' of radio sources — which we discuss in 
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We are now in a position to identify the timescale 



16 




2xl0 7 4xl0 7 



6xl0 7 8xl0 7 



lxlO L 



Age (years) 



Fig. 10. — Thick solid line: the trajectory of a point 
moving out from the centre of the AGN at which the 
break frequency has lowered to 1 GHz; thin solid line: 
the trajectory of one hotspot advancing away from the 
AGN; dashed line: the trajectory of the mid-point of 
the arm of a source (see text for more details). Upper 
plot: trajectories for a source whose jets each carry 
total kinetic energy Q Q = 10 40 W and lower plot: with 
Q = 5 x 10 37 W; the value of cj> was taken to be 5. 

8.2. Environments 

An important probe of the densities of the environ- 
ments of radio sources comes from studies of the de- 
polarisation of the polarised intensity from their syn- 
chrotron emitting lobes. Such a study by Garrington 
& Conway (1991) inferred the densities of particles re- 
sponsible for the depolarisation to be consistent with 



the notion that radio galaxies are in poor group envi- 
ronments, with a typical value of the density given by 
p = 1.67 x 10 _23 kgm~ 3 . We infer from Mulchaey 
& Zabludoff (1998) for density profiles of such groups 
which are parameterised as follows: 



p(r) = p 



-0 



(5) 



that suitable values of (3 and a Q are 1.5 and 10 kpc 
respectively (see also Willott et al. 1998b). 

Many authors have used various observational data 
as evidence for systematic changes in radio-source en- 
vironments with redshift. Examples include: (i) the 
X-ray detection of a number of high-z 3C sources 
interpreted as evidence that a larger proportion of 
these lie in rich (X-ray emitting) clusters (Crawford 
& Fabian 1996); (ii) a seeming increase in the optical 
richness of group/cluster environments with redshift 
(Hill & Lilly 1991, Yates, Miller & Peacock 1989); (in) 
the identification of 3C host galaxies at z ~ 1 with 
the brightest cluster galaxies (Best et al. 1998, Eales 
et al. 1997, Roche et al. 1998) and (iv) an increase 
in the radio depolarisation of radio sources with ei- 
ther increasing redshift or luminosity (Garrington & 
Conway 1991). 

We have attempted, and failed, to synthesise these 
observations into a consensus picture for the evolu- 
tion of radio-source environments. For example, the 
interpretation of X-ray detections associated with dis- 
tant powerful radio sources as emission from hot clus- 
ter gas by Crawford & Fabian (1996) is at odds with 
the interpretation of Brunetti et al. (1997) that in- 
verse Compton scattering from either the AGN radi- 
ation field or from the cosmic microwave background 
is responsible for these X-rays. While the finding of 
Garrington & Conway (1991) — that more distant 
sources show higher depolarisation — might appear 
to suggest that environments of more distant sources 
are more dense, it should be borne in mind that the 
more distant sources in their study, being shorter than 
lower- z sources, inform on the depolarisation of lobe 
emission at spatial points much closer in the AGN, 
where inhomogenities due to the host galaxy and the 
possible recent merger products which led to the for- 
mation of the AGN will increase the measured de- 
polarisation. In fact the interpretation of this very 
measurement is exacerbated by the "youth redshift" 
degeneracy, a point to which we return in § |9.5| . 

In the absence of any clear consensus on this ques- 
tion we have preferred to adopt the simple assump- 



17 



tion of no systematic change of radio-source environ- 
ment with redshift and to see whether, when imple- 
mented in our model for radio sources, this necessar- 
ily presents contradictions with observed properties 
of radio-source samples. One might expect informa- 
tion on radio-source environments to improve rapidly 
with the advent of facilities such as AXAF. 

8.3. Source sizes 

8.3.1. Source expansion 

The first essential step in modelling the time evolu- 
tion of the properties of an FR II radio source is to de- 
termine how the size of the source depends on its age 
t. Falle (1991) and Komissarov & Falle (1998) have 
shown that this dynamical problem can be solved via 
dimensional analysis, and since the mass involved is 
dominated by the ambient material swept-up within 
the bow-shock of the radio source the only character- 
istic quantities required are the kinetic power trans- 
ported by one jet Q a , and po x <Zq. This gives: 



D(t) = 2cia 



a 5 Po 



(6) 



where D(t) is the total length of a source, i.e. the 
hotspot to hotspot separation; see also Kaiser & 
Alexander (1997). Assumptions about the radio 
source environment are therefore necessary (as sum- 
marised in §|8.2|), but adopting these assumptions it 
is possible to constrain c\ in three ways: (i) for 
those powerful FRII sources with measured expan- 
sion speeds (v = 0.5dD/dt w 0.03c; Scheuer 1995, 
see §8.3.2) an upper limit may be obtained, similarly 
from (ii) the distribution of radio sources in the P—D 
plane of simulated samples (see Figures [l?] and |l^) 
and (iii) obtaining the axial ratios consistent with 
those observed. We take a value of 1.8 for the dimen- 
sionless constant c\ . Comparison of the observational 
and simulated data presented later in this paper (Fig- 
ures |l| k for the 3C data and Figures || & || for 
the 7C data) will be used to show that c\ must be 
similar to this value (see also Willott et al. 1998b). 
Falle's (1991) analysis also shows that the bow-shock 
must expand in a homologous (self-similar) manner. 

We wish to re-emphasize the distinction between 
the hotspot and the head region of the expanding ra- 
dio source (see §8.1.3 and Figure ft). The hotspot, also 



sometimes referred to as the 'working surface', may 
well move around the head of the lobe as in the "Den- 
tist's Drill" model of Scheuer (1982) so that it is very 



dangerous to use the (instantaneous) properties of the 
compact hotspot to determine how the (much larger) 
head region of the source expands, e.g. by equating 
the pressure of the hotspot to the 'ram pressure' of 
the external medium [e.g. Readhead et al. (1996)]. It 
is better to consider the expansion of the whole head 
region, and its associated bow shock. Similar points 
have been made by Begelman & Cioffi (1989). 

From the jump conditions at a strong shock, (see 
e.g. Landau & Lifshitz 1959) the pressure across the 
head (phcad) is related to the pressure in the external 
medium (p C xt) (whose adiabatic index is T K = 5/3) in 
terms of the Mach number of the bow shock (Mb) by: 
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In terms of the dimensional argument of Falle (1991) 
this may be expressed as: 
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where D is given by equation ^| and p cx t is the density 
of the external medium (identified with p(r) and given 
by equation ||). This, together with equation 0, gives 
an expression for the pressure in the head at time t: 
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Note that this is identical to equation 12 of Kaiser 
& Alexander (1997), but that we are making a clear 
association between environmental ram pressure and 
the average value of internal pressure across the entire 
head of the source, and not just the compact hotspot. 
It is significant in this regard that high resolution ra- 
dio imaging (Leahy et al. 1997, Hardcastle et al. 1997 
and Black et al. 1992) suggests that, ignoring the pos- 
sibly ephemeral hotspot emission, there is very little 
pressure variation across the head region of a given 
FRII radio source. It is also significant that the iden- 
tification of the characteristic pressure with the whole 
head of a source, rather than the hotspot, avoids the 
'lobe drag' problem posed, for example, by Williams 
(1991). 

We believe that the material backflowing from the 
hotspot, via the head, inflates a cocoon which we 
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identify with the observed radio lobe [an idea which 
dates back to Longair, Ryle and Scheuer (1973)] and 
which we assume confines the jet. Leahy & Perley 
(1991) made images at 1.4 GHz of 23 extragalactic 
radio sources with the VLA. In these images they de- 
tected the outer perimeter of the radio emission in 
most sources; at these outer edges the emissivity does 
not fade slowly into the noise, rather it suggests lobes 
which are over-pressured with respect to their local 



environment and will therefore expand outwards as 



radio sources end in shocks which mean that the ad- 
vance speed must exceed the ambient sound speed. 
The sound speed in an ionized gas at temperature T 
K is given by w 5 x 10~ 7 x c\fT where c is the speed 
of light; for T = 10 7 K this gives an advance speed of 
0.0015c. 

Upper limits to the advance speed come from con- 
sideration of the asymmetries in arm-lengths of radio 
sources. This kind of analysis was first done by Lon- 



thc source ages, and indeed as it expands longitudi- 
nally. 

Simulations show that jets with a low density ra- 
tio and high Mach numbers must produce backflow 
(Norman et al. 1982). Observational evidence for 
backflow, and high backflow speeds, comes from the 
spectral ageing analysis of Liu, Pooley & Riley (1992) 
which found lobe-speeds (i.e. the anti-vector sum of 
the oppositely directed head advance speed and the 
backflow speed) in some sources to exceed 0.2c. When 
this is taken together with Scheuer's (1995) analysis 
of arm-length asymmetries arising from light-travel 
time effects which limit the head advance speed to 
substantially below this, perhaps as low as 0.03c (see 
discussion of his method in § |B.3.2 ), the lobe-speeds 
which Liu et al. derive strongly point towards high 
backflow speeds out of the hotspot. This, in turn, is 
strong evidence for a pressure gradient from the head 
of the source back towards the central nucleus. 

If one considers a slice in low-frequency surface 
brightness through a radio lobe, along the jet axis, a 
gradual fall-off in surface brightness is seen from the 
head of the source back towards the core [see e.g. the 
slices shown in Leahy et al. (1989)]. Over a given vol- 
ume, pressure is related to surface-brightness under 
equipartition conditions by a simple power-law scal- 
ing: p oc S 4 / 7 . Thus we may infer a gradual decline 
in pressure from the head towards the core and that 
the head pressure governs the lobe pressure. 

The ratio of the pressure in the head to that in the 
lobe is not a well-defined quantity, but in the interests 
of simplicity in the modelling which we now describe, 
we will take for sources of all jet-powers a constant 
ratio of six for these pressures. 

8.3.2. Source advance 

A lower limit to the speed at which heads advance 
away from their AGN, and hence the rate at which 
radio sources grow, comes from the fact that FRII 



gair 



Riley (1979); a subsequent and similar study 
These studies assumed 



was done by Best et al. (1995) 



that all arm-length asymmetries arose from geomet- 
rical projection and light-travel time effects. No ac- 
count was taken of effects arising from environmental 
asymmetries. The former study obtained an upper 
limit to the head advance speed of 0.2c while the lat- 
ter study, which also took into account the effects of 
possible misalignment in the directions of ejection of 
the jets, obtained an upper limit of ~ 0.4c. These 
values should very much be regarded as upper limits; 
tighter constraints were obtained by |Scheucr (1995) 



who considered a sample of quasars each of which 
exhibited a considerable asymmetry in the presence 
of jets. The side of the brighter jet he ascribed as 
the nearer to us, attributing the jet brightness asym- 
metry to Doppler boosting. Arm-length asymmetries 
in these objects could then be more confidently at- 
tributed to light-travel time effects since there was 
a certain marker as to which was the nearer side to 
Earth. Scheuer's upper limit was somewhat sample 
dependent hinting at a non-constant hotspot advance 
speed, but gave a value of 0.03c for two samples and 
for all of his samples taken together a certain upper 
limit of 0.15c. This method of considering asymme- 
tries, gives an 'instantaneous' measure^] of the speed 
at which the heads were advancing out from the ac- 
tive nucleus at the time when the light we observe 
was emitted from the radio source. The term 'in- 
stantaneous' is used somewhat liberally here, to mean 
'within the light-travel time from the further hotspot 
to the nearer' as distinct from the speed which would 
be derived by averaging over the life-time of the source 
by dividing the (projected) length of the radio source 



5 The first direct measurements of source expansion were made 
by Owsianik & Conway (1998) from multi-epoch global VLBI 
observations; they measured an advance speed of the CSO 
0710+439 to be > 0.24c. In such a small source (its linear 
size is 101.5 pc) it is likely that the expansion is strongly 
determined by the local density profile within the core rariiuj 
of its host galaxy and so the environmental assumptions of f 
would be correspondingly poor. 
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by its total age. 

Speed limits in the same ball-park as Scheuer's 



smaller limits were obtained by Readhead et ah (1996) 



who investigated the empirical relationships between 
the external density and the size, pressure, advance 
speed and flux densities of compact symmetric objects 
by assuming the sources to be ram pressure confined. 
They argued that hotspot pressures adjusted to the 
ambient density so that head advance speeds could 
be constant and close to 0.02c (but see comments in 



.3.1) 



De Young (1997) has pointed out that the energy 
density at the end of the jet used to determine ad- 
vance speeds could be the result of cumulative ef- 
fects of the jet-environment interaction over the jet 
length. He suggests a preliminary approximation for 
how the advance speed depends on lossy jet propaga- 
tion, namely by considering a constant entrainmcnt 
rate per unit length of jet. He derives from this a 
roughly constant jet propagation rate. 

We chose to adopt the formalism for the linear size 
growth developed by Falle ( 1991) on dimensional ar- 
guments discussed in § B.3.1 . 

The jet-axes of the sources we simulate in our mod- 
els are deemed to be randomly oriented with respect 
to our line-of-sight. Each source is allocated an angle 
to the line of sight which is drawn from a distribu- 
tion uniform in (1 — cos 9). The projected linear size 
of each simulated object is derived by evaluating D, 
the absolute physical size of the radio source, given 
in equation |^ at the time when the source intercepts 
our light-cone, and multiplying this length by sin 8. 

In our model we deem that all radio sources will 
have their beams switched on for some maximum 
length of time (i m ax-agc, which has been taken as 
5 x 10 8 years). The time at which we 'make the obser- 
vation' (i.e., examine whether the source is above the 
flux-limit of our simulated survey) of a given source 
is randomly chosen from a distribution described in 
§|1| between years and i max _ ag o years. If £ max _ ago 
were chosen to be too short, there would be two erro- 
neous consequences: (i) if the true value of the coef- 
ficient of proportionality for the head advance (ci) is 
known then sources will never achieve projected linear 
sizes as large as those which are observed and (ii) a 
distribution of points would be generated on the P-D 
plane, with an abrupt maximum size across the whole 
range in power, strikingly different from that which is 
observed. This would arise because the sampling of 



the radio sources would occur at early stages in their 
lives when very few sources would have fallen below 
the sample flux- limit. If the sampling of the radio 
sources broadly occurs over a time when the 'natural 
selection' process (the catastrophic decrease in lumi- 
nosity due to synchrotron, adiabatic expansion and 
inverse Compton losses) governs the distribution of 
points on the P-D plane, then it much more closely 
resembles that arising from real data. 

If imax-agc is chosen to be too large, then provided 
the simulations are allowed to run through the birth of 
a sufficiently large number of radio sources there will 
be no consequence on the distribution of points on the 
P-D plane. This is the case because the sampling of 
radio sources during the period of their lives when 
they are sufficiently bright to be detected by a survey 
on Earth is still random throughout all of this period. 
That there may be a relatively long period of time 
when radio sources have fallen below a sample flux- 
limit during which our light-cone randomly samples 
them only has consequences for the number of sources 
which have to be 'observed' until the simulated survey 
is complete. 

8.4. Luminosity and spectral index develop- 
ment 

We took as our basic model for the luminosity- 
with-time profile of a given radio source the model 
developed by Kaiser fc Alexander (1997)| and Kaiser 



Dennett-Thorpe & Alexander (1997); this is the first 



model of extragalactic radio sources which directly 
incorporates the effects of inverse Compton losses on 
the luminosity of a source. This model assumes that 
source expansion is described by Falle's (1991) model, 
i.e. equation ^. 

8.4-1- The role of hotspots in governing the lobe lu- 
minosity 



As described in §3.7 each object in our model is 



assigned two jets which each transport a bulk kinetic 
power of Q , assumed to be constant throughout the 
lifetime of the source. 

Under the assumptions of 



3.1.2 



we 



can assume 

that the thrust of the jet is given by (Q a /c). The 
stagnation pressure p s in the post-jet shock material 
will then be approximately given by: 
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where Ah s is the area perpendicular to the jet over 
which the jet thrust acts. We equate this stagnation 
pressure with phs, the pressure in the hotspot itself. 
If we assume that the energetics of the hotspot are 
described by equipartition, then the total energy den- 
sity in the hotspot is due to two equal contributions, 
one from the particles and one from the magnetic field 
(this latter being equal to B^ a /2fi ). We may then ex- 
press the magnetic field (assuming it to be tangled) in 
the hotspot (-Bhs) m terms of the jet power as follows: 



c 3^ 



(11) 



since for a relativistic fluid, and also for a tangled 
magnetic field, pressure is one-third of the total en- 
ergy density. The break frequency in the hotspot vy^, 
is inferred from the usual expression for the break 
frequency in a synchrotron source: 



(9/4)c 7 B hs 
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where following Leahy (1991), the coefficient of pro- 
portionality c 7 is 1.12 x 10 3 nT 3 Myr 2 GHz thus the 
hotspot magnetic field Bhs and equivalent magnetic 
field due to the cosmic microwave background Bqmb 
are measured in nT, t s the 'synchrotron age' of the 
distribution is measured in Myr and the break fre- 
quency of the hotspot isbh is measured in GHz. In 
the context of our model, we assume that the syn- 
chrotron age of the radiating particles depends on how 
long the particles are exposed to the magnetic field of 
the hotspot before they escape into the lobe. Some 
effect of this type seems inescapable in any model 
for a hotspot [e.g. Heavens & Meisenheimer 1987 and 
Eales, Alexander & Duncan 1989]. The magnetic field 
has the dependence on Q given by equation [ll] giving 
a value for the break frequency which depends just on 
Qo, if we make the assumption that hotspot sizes are 
constant (we assume for all our sources that the ra- 
dius and depth of each hotspot are 2.5 kpc and 1 kpc 
respectively). 

8.4-2. Lobe luminosity 

Our formalism essentially follows the model of 
Kaiser et al. (1997) but with two important differ- 
ences: (i) instead of assuming a constant injection 
index [given by p — —d\ogn{^)/d\og^f\ of p = 2.14 
we assume that the injection index is governed by the 



breaks in the energy distribution of the particles in- 
jected into the lobe from the hotspot |] These are ul- 
timately determined by the maximum and minimum 
particle dwell-times in the hotspot together with the 
magnetic field strength in the hotspot, which as seen 
in equation [n] directly depends on Q Q itself, (ii) We 
also assume that adiabatic expansion losses out of the 
hotspot are governed by the hotspot pressure as de- 
rived in equation [l^, not by the pressure in the head 
of the source which governs the increase in arm-length 
of the source. 

Note that this process of adiabatic expansion of 
the particles out of the hotspot into the lobes does 
not move the energy distribution of the lobes out of 
equipartition, if the adiabatic indices of the particles 
and of the fields are both 4/3 (as is the case for a 
highly relativistic fluid and for the virtual photons of 
a magnetic field). 

A continuum of spectra may be regarded as in- 
jected to the lobe as the supply of particles in the 
hotspot is tapped for those populations of particles 
which have spent different dwell-times in the hotspot 
and thus have different break frequencies. Since the 
magnetic field in the hotspot is governed by the jet 
power Q Q , this ultimately governs the spectrum of 
particles injected into the lobes [this is the origin of 
the P-a correlation in our model]. 

We derive a break in the energy spectrum of syn- 
chrotron particles freshly injected into the lobe by 
consideration of the corresponding break in their fre- 
quency spectrum of synchrotron emission. Since spec- 
tra are unchanged in shape, though translated to 
lower amplitudes and lower frequencies, by adiabatic 
expansion (Schcucr & Williams 1968) — e.g. from 
the hotspot into the lobe — the Lorentz factors of 
the break frequencies of the particles injected into 
the lobes at some time ti n j may be derived. Scheuer 
& Williams showed that if the factor by which the 
linear scale of the emitting population has changed 
after undergoing adiabatic expansion is /, then the 
factor relating the break frequency after expansion to 
that before is /~ 4 . It is assumed that at the instant 
that an elemental population of particles is injected 
from the hotspot into the lobe, they undergo expan- 
sion which is approximately adiabatic. The break fre- 
quency of the synchrotron population as it is injected 
into the lobe is related to its break frequency in the 



6 We remind the reader of the distinction we draw be 
jection index and acceleration index as discussed in §5 



ween m- 
1. 



21 



hotspots by a factor of (phs/piobe[*inj]) -4 ^ 3r \ where 
(phs/piobe^mj]) 1 ^ 31 ^ is their linear scaling factor (r 
is the usual adiabatic index of 4/3) and where t; n j 
is the time at which a population of particles is in- 
jected into the lobe. This may be expressed in terms 
of the magnetic energy densities instead of the pres- 
sures i.e., the break frequency of the synchrotron spec- 
trum of the particles freshly injected into the lobes 
is (B^ s /(2^oUs[tinj])) _4/(3r) x i/ hh where h> hh was the 
break frequency in the hotspots and ws^mj] is the 
magnetic energy density in the lobes at this time ii n j 
(which is directly proportional to the pressure in the 
lobe and head given in equation ||) . 

The energy of such a particle, freshly injected into 
the lobe, which is emitting at the break frequency v^i 
is therefore given by 7biwi c 2 where 7bi is given by: 



log n("/) 



7bi[i 



in.) J 



^bi[*i] 



^bh 



2/j, u B [th 



(13) 

where u g is the non-relativistic gyro-frequency, given 
by eB[t inj }/(2nm c ). 

The exponent of the energy spectrum in the low 7 
regime (pi) is taken to be 2 (corresponding to a fre- 
quency spectral index a of 0.5) assuming the particle 
acceleration mechanism in the hotspots to be via the 
first-order Fermi mechanism (Bell 1978). We note 
however that it is quite likely that the value of the 
acceleration index might be lower than 0.5 (Heavens 

6 Drury 1988) although we do not implement values 
lower than this in our model. 

The injection spectrum is derived from the two 
break frequencies in the hotspot emission at 24, s and 
Vhf- The former is the point in the spectrum where 
the population of particles which has lingered for the 
longest time in the hotspot changes from a spectral 
index of a = 0.5 to 1 (or p = 2 to 3, since p = 2a+ 1). 
The latter is the point in the spectrum where the pop- 
ulation of particles which escapes most quickly from 
the magnetic field of the hotspot changes from a spec- 
tral index of a = 0.5 to 1 (or p = 2 to 3). 

We deem the maximum 7 (i.e. 7b s ) in the p = 2 
regime to be the Lorentz factor corresponding to the 
break frequency of the electrons which linger longest 
in the strong magnetic field of the hotspot before be- 
ing injected into the lobes. 

The exponent of the energy spectrum in the high 

7 regime (ph) is taken to be 3 in identifying a hotspot 
as an example of continuous injection (Carilli et al. 
1991). We deem the minimum 7 (which is 7b f) in 
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Fig. 11. — This schematic illustration shows the dif- 
ferent energy regimes within a hotspot. The low en- 
ergy spectrum is described by a region with index 
p = 2 (in light grey). The high energy spectrum is 
governed by a region with index p — 3 (in dark grey) . 
The middle region is described by a single power law 
with index p = p m (in white). The precise value 
of this index depends on how leaky the hotspot is 
deemed to be and the difference in dwell-times in the 
hotspot of the particles which leak out most quickly 
and those which leak out most slowly. 



this regime to be the Lorentz factor corresponding 
to the break frequency of the electrons which leak 
out most quickly from the strong magnetic field of 
the hotspot into the lobe (see Figure |ll|) . The shape 
of the spectrum between 7b s and 7b f is determined 
by the relative rates at which particles with different 
dwell-times in the hotspot leak out into the lobe. 

A refined model of the details of this mechanism 
might suggest that this region should have quite a 
curved spectrum, but in the interests of simplicity, 
we approximate this region by a straight line whose 
gradient can take a maximum value of 3 (a — 1) and 
a minimum value of 2 (a = 0.5). 
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The values of 7b s and 7b f are derived by considering the longest and shortest times for which particles reside 
under the influence of the magnetic field of the hotspot and deriving the appropriate break frequencies using 
equation [l^. These dwell-times were taken to be 10 5 years and 1 year respectively corresponding to, in the former 
case a backflow speed of 0.03c and, in the latter case, to particles which do not fully cross the hotspot. The break 
frequencies from these may be converted to Lorentz factors corresponding to breaks in the energy spectrum using 
equation [l3| The derivation of Ub{= ^[Tbs]) (the number density of particles with a given Lorentz factor 7b s per 
unit Lorentz factor, see Figure at a time U n j is obtained by equating the energy density at time ti n j with the 



integral over the whole energy range, from 7 n 
follows: 



to 7 m ax for each particle of energy 7 appropriately weighted as 



u c [t 



mjj 



7i h (7)7m e c 2 d"f ■ 
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(14) 



with the weighting factors 711(7), n m (j) and rih("f) from the three regimes in 7 respectively: 

-2 



7 

711(7) = n b[ 

\7bs 



n m (7) = n b 

\7bs 



if 7 < 7b s, 



if 7bs < 7 < 7bf, 



(15) 



«h(7) = f n b if 7 > 7b f 

V7bf/ 

where / = (7bs/7bf) Pm , as shown in Figure [LT[ 

Thus we find that the normalisation factor is given by: 
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(16) 



The energy density in the lobe at time t, u c [t], is directly proportional to the pressure in the lobe and the head at 
time t, given in equation ^, with coefficients of proportionality which are of order unity (see Kaiser et al. 1997). 

We take 7 m ; n to be 1 and j max , the maximum value of the Lorentz factor in the hotspots, to be 10 14 since 
emission in the intense magnetic fields of the hotspots is observed to extend to frequencies of several tens of GHz. 

We now consider the thermodynamics of the transition of the relativistic plasma from the jet, to the hotspot 
(meaning the post-jet shock region), and then through the head to its eventual home in the lobe itself. Following 
Kaiser et al. (1997) we equate the final energy state of an elemental volume 8V\ of plasma in the lobe with the 
energy it had when the jet delivered it to the hotspot (a dominant fraction of which is assumed to be in the form 
of bulk kinetic energy) minus the work done in adiabatically expanding from the hotspot to the lobe. We assume 
that the jet-power Q a is completely converted into thermal energy at the hotspot. This seems reasonable given 
that the flow speeds into hotspots are close to c (Wardle & Aaron 1997) while those out of hotspots can be up to 
~ 0.2c and are oppositely directed (Liu et al. 1992). [The energy lost to radiation mechanisms is neglected and can 
be shown to be many orders of magnitude smaller than the other quantities involved.] Hence, 
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= QoSt 
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(17) 



where 8V\ is the volume element injected into the lobe, <5Vhs is the volume which this lobe clement occupied when it 
was in the hotspot, <5ti n j is the time interval at tj n j over which this element is injected into the lobe, p^s and p\ are 
the pressure in the hotspot and lobe respectively and T\ and Ths are the adiabatic indices of the lobe and hotspot 
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respectively each taken to be 4/3. Since the expansion out of the hotspot is approximated to be adiabatic, we may 
write that: 



Phs SVhs ( p h ; 



pi SV\ \ p\ 



(18) 



and using equation ml we may thus write the contribution to the volume of the lobe from the element injected in 
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After multiplying by a term to allow for the further (adiabatic) expansion of the element as the source evolves 
(namely [i/^mj] ai , where a\ is a constant which relates the expansion of volume V after time t as V oc t ai ), the 
integral of equation [l^ over t- m j may be identified as the volume V\{t) now occupied at time t which has been 
injected into the lobe since time t m i n , and hence as the volume term in equation @; this is given by: 



(Tj - 1)Q ( Qo 



dt iai . (20) 



Note that the net result of this expansion is that only a fraction \i of the thermal energy deposited in the hotspots 
by the jet appears as stored energy in the lobes, the remainder being lost as work done against the shocked ambient 
medium in the bow-shock. The magnitude of /i can be estimated from equation [l8] which is is essentially the 
reciprocal of the ratio of the internal energy of a given element in the lobe to that which it had in the hotspot. For 
Ti = Ths = 4/3, and typical values of Phs/piobc ~ 10 - 100, /i varies between 0.6 - 0.3. 

Figure [l2] indicates the six possible ways in which the slope of the energy distribution of the relativistic particles 
injected into the lobe (i.e. the gradient at 7i n j) may change with time, according to the relative value of 7i n j 
compared to 7b s and 7b f . 

We now consider the relation between the Lorentz factor of a particle injected into the lobe which will later 
contribute to the emission we observe at i bsi an d its Lorentz factor at the time when the particle is injected. 
At £ bs (the point when the light we observe leaves the object, i.e. when the object intercepts our light-cone) the 
Lorentz factor of the particles injected at that point in time, must be given by equation |3| for a given emitted 
frequency. The particles which contribute to the emission leaving the object at t ODS but which were injected into 
the lobe at a time At prior to t ODS must have had a Lorentz factor higher when they were injected than they have 
at i bs, to compensate for the energy depletion during At. The maximum value of At gives i m i n , the time before 
which no injected particles contribute to < bs- As in Kaiser et al. (1997) t mm is the time before which no particles 
which give rise to the emission observed (at the chosen frequency) at t oos can be injected because the losses they 
suffer, given their Lorentz factors, would be too catastrophic. Thus for the ensemble of particles contributing to 
the radiation at a given emitted frequency at £ bs, those with the largest Lorentz factors at injection are injected 
at t m i n and those with the lowest are those actually injected at t oos . 

We now consider the relation between the Lorentz factors injected at time i m i n and t oos and the corresponding 
breaks in the injected energy spectrum (depicted in Fig. |ll|). If the Lorentz factor of the injected electrons is 
always above these two or always below or always between them then the calculation is straightforward involving 
just one integral with the weighting function in each case being just either the bottom or the top or the middle form 
respectively of the full form of the weighting function listed in equation [l5]. If, during the time between i m ; n and 
tabs the gamma-break should overlap with two distinct regimes in the spectrum, the integral must be performed in 
two steps. If during the time between t m i n and t bs the gamma-break should overlap with all three distinct regimes 
in the spectrum, the integral must be performed in three steps (see Figure [l2|), thus: 
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Fig. 12. — Schematic depiction of the six possible ways in which the Lorentz factor at injection of those particles 
which ultimately contribute to the radiation we observe (emitted at i bs) may relate to the Lorentz factors at the 
breaks in the injected spectrum. As in Figure [ll], the light grey regions are where p = 2, the dark grey regions are 
where p = 3 and the white regions are where p = p m . 
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where the volume element in equation |2j has been replaced by the integrand on the right-hand side of equation [20|, 
and where itB^inj] represents the magnetic energy density in the lobe at the time a given element of particles is 
injected. The pressure in the lobe p\ is taken to be a constant factor of six lower than the pressure in the head 
Phead given by equation ||. We use Kaiser et al.'s (1997) expression for 7[ii n j] in their equation 10. 
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8.5. Predictions for the evolution of individ- 
ual sources 

8.5.1. Sizes and energy budgets 

The volume of the radio source at the time of 'ob- 
servation' may be found by integration of the volume 
elements which have been injected over the life of the 
radio source from a very early time in its life (e.g. 
fmin = 1 x 10 years), up to t a h s , i.e. by performing 
the integral of equation |20j, to allow for the expansion 
of these elements between the time they were injected 
and the time the 'observation' is made. 

The energy transported, E%, by the jets in a radio 
source, after they have been switched on for a time t, 
is given simply by 



E t — 2Q D t. 



(22) 



The energy, E s , which is stored in the lobes according 
to the usual minimum energy arguments depends on 
the luminosity and the volume of the lobes as follows: 



R 



6/7 



(23) 



where x is a constant throughout the life of a source 
which depends on assumptions about the low-energy 
cut-off of particles and their filling factor in the ra- 
dio lobes (see Miley 1980), P151 is the luminosity at 
rest-frame 151 MHz, D is the total linear size of the 
radio source and where Rr is the axial ratio of the 
radio source, defined as the total length divided by 
the width, as one would measure on a radio map. 

In our model, throughout a source's life, it does 
not store energy at a rate which is proportional to its 
age; this arises because the rate at which work is done 
by the expanding source itself depends on the age of 
the source. This follows from equation |l8| and the 
fact that hotspot pressure does not scale with time 
in the same way as lobe pressure (unlike self-similar 
models). The prediction of axial ratios by our model 
comes from finding (the square root of 7r/4x) the 
quotient of the cube of the linear size derived from 
equation |g| and the volume obtained from the inte- 
gration of equation ^Ol We find that throughout the 
lifetime of an individual source its axial ratio steadily 
increases thus, for an individual source its expansion 
is not self-similar. For sources with Q = 10 38 W, 
we find that the axial ratios increase from 2 to 8, 
while for sources with Q = 10 40 W, the axial ratios 
increase from 4 to 13. Our model thus reproduces 



axial ratios which are larger in sources with higher 
jet-powers, consistent with the findings of Leahy & 
Williams (1984) and Leahy et al. (1989) who found 
axial ratios between 4 and 12 (using the definition 
that the axial ratio is the total projected length of 
the source divided by its width). We note that the 
closer the jet- axis is to the line-of-sight, the smaller 
the measured axial ratio of a source will be because of 
projection effects. Both Leahy & Williams and Leahy 
et al. found evidence of a trend of increasing axial ra- 
tio with increasing source luminosity and commented 
that the scatter was greater than that expected from 
projection effects alone. Leahy & Williams' study was 
for a subset of the 3C sample (see §3T) with angular 
sizes greater than 45". With the elimination of these 
smaller (hence, given the sample, younger) sources it 
is then unsurprising in the context of our model that 
they found their more powerful sources to have higher 
axial ratios. 

8.5.2. P-D tracks 

Figure [l3] shows the tracks across the P-D plane of 
sources according to our model described above with 
various values of jet-power and redshift and compares 
them with the predictions of the model of Kaiser et al. 
(1997). Note that in each case the tracks may be ap- 
proximately described by two power-laws, either side 
of the knee. In each case on the left-hand side of the 
knee the gradient of our tracks are steeper than those 
of Kaiser et al. This arises because in our model the 
synchrotron emitting particles adiabatically expand 
out of a constant pressure hotspot into a lobe whose 
pressure falls throughout the source lifetime, thus the 
losses due to adiabatic expansion increase with time. 
In the model of Kaiser et al. (1997), the expansion of 
the synchrotron emitting population is from a head 
into a lobe whose pressures are in a constant ratio 
throughout the lifetime of the source. Thus the adia- 
batic losses in this latter model do not increase with 
time. The increasing adiabatic expansion losses with 
time in our model ensure that the tracks in the P-D 
plane can be very steep even in an environment whose 
(3 exponent is not particularly large (e.g. 1.5 here as 
described in § 3.2) . This point plays a crucial role in 
determining the ease with which (particularly large, 
powerful) sources fall through the survey flux- limit. 

The 'knee' between the two regimes of a P-D track 
occurs earlier in the lives of sources at higher red- 
shifts. Below the knee the magnetic field in the lobes 
-Biobc is stronger than the equivalent field due to the 
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Fig. 13. — We overlay on this plot in solid lines 
the tracks through the P-D plane of three sources 
(upper to lower) with Q a = 1.3 x 10 40 W at z = 2, 
Q = 1.3 x 10 39 W at z = 0.5, Q Q = 1.3 x 10 38 W at 
z = 0.2. We overlay in dotted lines the tracks from 
the model of Kaiser et al. (1997) for sources with the 
same jet-powers and redshifts. The larger symbols in 
this plot show the rest-frame luminosities at 151 MHz 
of members of the 3C complete sample (§|3.l|) against 
their projected linear sizes. Since this is a logarithmic 
plot, the differences in the distributions between the 
projected and the true linear sizes of sources will be 
slight. Small dots lie on the tracks to indicate the 
times 1, 10, 20. . . 100, 200 Myr through the source's 
life. All tracks in this figure wer e ge nerated assuming 
an environment as described in §8.2, i.e. with [3 = 1.5 
and with c\ = 1.8. 



cosmic microwave background Bcmb ■ Spectral losses 
for a particle of Lorentz factor 7 go as 7 2 -Biobe> but 
for a given 7 this rate of energy loss decreases with 
time since i?i bo is decreasing with time. Through- 
out the lifetime of the source, Biobc will be falling — 
thus for radiation at a fixed frequency which we might 
observe in our survey, the Lorentz factor 7 will be pro- 
portional to the reciprocal of the root of the lobe mag- 
netic field (see equation |J). Thus as the source ages, 
increasingly energetic particles are responsible for the 
emission observed at some fixed frequency. Once the 
magnetic field of the lobe has dropped below -Bomb, 
the inverse Compton losses will dominate the rate at 
which the high 7 particles are depleted, in this sec- 



ond regime in the source's life. This magnetic field is 
constant and not decreasing as -Biobo is, so the spec- 
tral losses thereafter remain at a fixed, not a decreas- 
ing, rate. However, the lobe magnetic field i?i bo con- 
tinues to fall requiring yet higher Lorentz factors for 
emission at a given frequency, and the high Lorentz 
factor population is being depleted at a steady rate 
due to inverse Compton scattering from the cosmic 
microwave background. 

Figure [li] shows tracks across the P-D plane de- 
rived from our model for various values of jet-power 
and redshift. The solid lines of the tracks become dot- 
ted when the luminosity and redshift of that source 
are such that they fall below the survey flux-density 
limit of 0.5 Jy. Comparison of the time markers 
along the different tracks show that the tracks of 
the more powerful, higher-redshift sources become 
steeper, this together with the consequences of the 
flux-limit, means that these sources fall below a sur- 
vey flux-limit at an earlier stage in their lives than 
do sources at lower power and redshift. Here then, is 
an indication o f the origins of the effect we discuss in 
more detail in §|9.5| , the "youth-redshift degeneracy" 
which is responsible for linear size evolution, on the 
assumption that radio sources are intercepted by our 
light-cone in a manner which is essentially random. 
We examine in detail the correct way to calculate the 
interception of the source population with our light- 
cone in §8.6 and then see clearly how these two effects 
arise. 

8.6. Redshift distribution 

In this section we consider how to model the num- 
ber of radio sources born at a sequence of cosmic 
times, and their subsequent interception with our 
light-cone. Successive populations of radio sources 
are hatched every 10 6 years, and are evolved, starting 
at some initial high redshift [z = 10 is used). 

We assume that the number of sources born at a 
given cosmic time (t), per unit cosmic time, per unit 
co-moving volume element (defined as having a proper 
volume of one cubic Gpc at redshift zero) is given by 
p(t) dt = p(z) dz which may be parameterised by: 



p(z)dz oc exp ^ ( Z ^* 1 



dz. 



(24) 



We further assume that at any given cosmic time, 
this description holds at any spatial point in the 
Universe, since this is presumed to be homogeneous 
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Fig. 14. — We overlay on this plot the tracks through 
the P—D plane of five sources (upper to lower) with 
Q = 5 x 10 38 W at z = 3, Q Q = 2 x 10 38 W at 
z = 1.5, Q = 2 x 10 38 W at z = 2.5, Q a = 1 x 10 38 
W at z = 0.8, and Q = 5 x 10 37 W at z = 0.5. The 
left-most light-grey bar indicates that for an object 
to be brighter than the flux-limit of the 7C sample 
(0.5 Jy), if it is at z > 2 then its luminosity must 
be in the range indicated by the bar; the correspond- 
ing bars going to lower luminosities are for the lower 
redshifts 1.5, 0.8 and 0.5 respectively. The solid lines 
of the tracks become dotted when the luminosity and 
redshift of that source are such that they fall below 
the survey flux- limit of 0.5 Jy. The larger symbols in 
this plot show the rest-frame luminosities at 151 MHz 
of members of the 7C complete sample (§|3.3|) against 
their projected linear sizes. Small dots lie on the 
tracks to indicate the times 1, 10, 20. . . 100, 200 Myr 
through the source's life; these become cross-symbols 
when an object falls below the flux-limit. It can be 
seen that because the tracks of the more powerful, 
higher-redshift sources become steeper, and because 
of the consequences of the flux- limit, these sources fall 
below a survey flux-limit at an earlier stage in their 
lives than do sources at lower power and redshift. The 
same environment and value of c\ are assumed as for 
Figure 13. 



and isotropic. We make the approximation that all 
sources are precisely following the Hubble flow. 

We wish to consider the total number of sources 
born at a given cosmic time which will be intercepted 



by our light cone, within a certain interval in cosmic 
time from the time when they were born. [This inter- 
val is chosen to be the maximum length of time for 
which the beams in each radio source are switched on. 
We assume that in each radio source this is 5 x 10 8 
years (see justification of this number in § |S.3.2 ), pre- 
suming no dependence on redshift or on intrinsic jet- 
power.] If a source born at a given point in cosmic 
time is to be intercepted by our light-cone at some 
point between the time when it was born and the time 
when its beam is switched off this has implications for 
the distance it must be from us at the point when it 
is born. If the source is intercepted by our light-cone 
at the time when it is born, then its distance from us 
is simply given by the distance along our light-cone to 
this redshift (which is hereafter referred to as Zbirth)- 
If the source is intercepted by our light-cone at the 
moment when the beam is switched off (i.e. 5 x 10 s 
years after the source was born) then its distance from 
us at this later point in cosmic time must be the dis- 
tance along the light-cone to this lower redshift, but 
this must correspond to a specific distance from us 
at the earlier cosmic time when the source was born. 
With this information a thick radial shell of a specific 
proper volume whose maximum and minimum radii 
are defined by the maximum distance and the mini- 
mum distance a source may be from us if it is born 
at a given ^birth and if it is intercepted by our light- 
cone within a cosmic time of 5 x 10 8 years from Zbirth- 
This region, for a given Zbirth, is hereafter called the 
'relevant co- moving volume element'. The relevant 
co-moving volume element is the maximum size, for 
a given solid angle, of a box which contains all the 
sources which must be intercepted within an interval 
of 5 x 10 8 years of being born at the redshift Zbirth 
(and no sources born at this redshift which are not 
intercepted by our light-cone within this interval). 
The form of the metric which we use in this analysis 

is: 



ds 2 = c 2 dt 2 - 



R 2 (t) 



(1 + i£x 2 ) 2 



[dr 2 +r 2 (d0 2 



(25) 

where R(t) is the scale factor of the Universe at cosmic 
time t, r is a radial co-moving co-ordinate, 9 and 4> are 
the polar and azimuthal co-ordinates; k refers to the 
curvature of the Universe (k = for a spatially flat 
Universe, k = 1 for an open Universe and k = — 1 for 
a closed Universe). Following standard practice, we 
introduce an alternative radial co-moving co-ordinate 
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X defined in terms of r by: 



X = r if k = 

sinx = r/[l + r 2 /4] if k = 1 
sinhx = r/[l — r 2 /4] if k = — 1 



(26) 



writing the space-time metric as 

ds 2 = c 2 dt 2 - dl 2 , 
we have for k = 0: 

dZ 2 = i? 2 (t)[dx 2 + X 2 {d0 2 + sin 2 6># 2 )] 
while for k = 1, we write: 

rf/ 2 = i? 2 (£)[dx 2 + sinx 2 ^ 2 + sin 2 6»d0 2 ) 
and for k = — 1, we write: 

dl 2 = R 2 [t) [dx 2 + sinh x 2 (^ 2 + sin 2 (9d</> 2 )] . (30) 



(27) 
(28) 
(29) 



In the three cases the proper volume A V of a spher- 
ical shell, whose inner and outer edges are at xi and 
Xi respectively, (a co-moving volume element) is given 
by: 



AV 



2irR 3 (t) 
2nR 3 (t) 
2nR 3 (t) 



3* 



X2 



Xi 



X- 5sin2x 



for k = 



for k 



i sinh 2x - X 



Xl 
Xl 



for fc 



Xl 



(31) 

Along a light-ray, the interval ds — 0, and along a 
radial light-ray 

ds 2 = c 2 dt 2 - R 2 (t)d X 2 (32) 

is independent of the value of k so we may write 

1 



X 



R(t) 



dt. 



(33) 



The limits of integration, xi an d Xii are respec- 
tively the radial co-ordinates to the outer and inner 
surfaces of the thick radial shell and are given by: 



Xi 



(no 



1 



dt. 



and 



Xi 



tbirth+5xl0 8 yr 



R(t) 



dt. 



(34) 
(35) 



where tbirth is the cosmic time when this batch of 
sources is born. 

In order to convert proper volumes measured in 
cubic Gpc at the time of birth into units which will 
be cubic Gpc at redshift zero, the proper volumes are 
multiplied by (1 + Zbirth) 3 - 
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Fig. 15. — This figure shows the volume of the 'rel- 
evant co-moving element' (see text for definition; the 
units are those which will be cubic Gpc at redshift 
zero) at a given redshift, in each of the three different 
cosmological models we consider in this paper. 

Once the volume of this co-moving element is cal- 
culated in a chosen cosmological model, it is then sim- 
ple to calculate the number of sources born [per solid 
angle born in the chosen interval in cosmic time] at 
a given Zbh-th which will ultimately be intercepted by 
our light-cone, by multiplying this by p{z) given in 
equation ^J. Precisely when each source intercepts 
our light-cone is determined by its radial distance 
from us at the cosmic time corresponding to Zbh-th 
within the radial shell described above. Since we 
assume spatial homogeneity we deem these sources 
to be randomly distributed within the volume of the 
box. Thus the length of time which elapses after a 
given source is born before it intercepts the light- 
cone, (called t a ge) is randomly drawn between and 
5 x 10 s years, but not from a uniform distribution. 
The uniformity of space requires that within equal 
(infinitesimal) volumes within the thick radial shell 
there should be equal numbers of objects born at 
Zbirth- Sources are randomly allocated to equal vol- 
ume shells, then their radial co-ordinate (x) is de- 
rived by inverting the equation [3l] appropriate for the 
value of k chosen, and solving this using an iterative 
method. We see that this process means that a source 
is therefore more likely to have a larger radial co- 
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ordinate within a given spherical shell than a smaller 
one, since the sources are uniformly distributed in vol- 
ume rather than in their radial co-ordinate, and hence 
there is a mild bias arising from geometry for a source 
to be seen at an earlier rather than a later stage in 
its life. Once the appropriate radial co-ordinate is 
drawn, the cosmic time (t bs) at which the source is 
intercepted by our light-cone may be determined from 



X 



1 



R(t) 



dt. 



(36) 



By subtracting the cosmic time when the source is 
intercepted by the light-cone from the cosmic time 
when the source was born, the time for which a source 
should evolve i.e. t agc , according to its initial param- 
eters, is known. Thus, 



tr 



^birth- 



(37) 



Therefore, at a given point in cosmic time, we 
know the exact number of sources born which will 
intercept the light-cone thereafter within a maximum 
time-period of 5x 10 s years, and the length of time for 
which each evolves until the light ultimately observed 
by us is emitted. The redshift of such an object is then 
derived by considering the cosmic time at which the 
source was born, plus this time for which the source 
lived (t age ) until the light which we ultimately observe 
on the Earth, left the source. 

8.7. Parameters and results 

A jet-power is allocated to each jet in a simulated 
object which has been drawn randomly from a distri- 
bution of the form: 



p(Qo)dQ Q = 



Qo if Qmin ^ Qo ^ Qmaxi 

if Qo > Qmax or Q Q < Q mi 

(38) 

where x is a positive number. Each source is then 
allowed to evolve in lu min osity (according to the pre- 
scription discussed in § 3T ) at a (rest-frame) frequency 
( 1 + z) x the observing frequency for the length of time 
< ag0 . If the luminosity of the source at this frequency 
is such that its flux is above the deemed survey flux- 
limit (e.g., 0.5 Jy for a simulated 7C sample) given its 
redshift, it is regarded as being detected by the sur- 
vey. Sources which do not meet the survey flux-limit 
are thereafter ignored. Note that we do not simulate 
the effect of surface brightness, as distinct from flux 
density, in these simulated surveys. However, as Ri- 
ley (1989) discusses the numbers of sources omitted 



from real surveys due to surface brightness effects in 
large angular size sources is small indeed. 

For those sources which do lie within the sample 
flux-limits, the following calculations are performed: 
the luminosities at two other frequencies (includ- 
ing rest-frame 151 MHz), for the same point in the 
source's lifetime, are derived as is the projected lin- 
ear size (according to the prescription discussed in 
§8.3.2). From three (luminosity, frequency) points 
now derived for each source in the sample, a quadratic 
expression may be fitted from which the spectral in- 
dex at rest-frame 151 MHz is derived, together with 
the a\ and ai coefficients of spectral shape discussed 
in Section 



4.2 



Figures 



It a nd |l7| show real and simulated 3C data. 
Figures |18 and |19| show real and simulated 7C data. 
In these simulations the 'sky-areas' (i.e. the solid an- 
gle covered by each survey) differ by a factor of 430, 
which is of the same order as the factor by which the 
true sky areas covered by 3C and 7C differ (4.2 sr 
versus 0.013 sr). For the simulations presented here 
the values of the Q -related quantities included in 
equation |38| are x = 2.6, Q max = 5 x 10 42 W and 
Qmin = 5 x 10 37 W. The values of the z-related quan- 
tities are z\ = 2.2 and Az% = 0.6. 

8.8. Discussion of simulations 

A comparison of the simulated and real 3C data 
suggests that to first order there is agreement in the 
basic quantities. Sources are detected out to ap- 
proximately the same redshift, the maximum size of 
sources detected is approximately correct, the depen- 
dence between linear size and redshift shows encour- 
aging correspondence and spectral indices increase 
with luminosity rather than with redshift. These cor- 
respondences also seem to be broadly present in the 
7C datasets though there is of course considerably less 
scatter seen in the spectral index distributions; note 
that we have not included in our model the possibil- 
ity that acceleration indices within the hotspot might 
be less than 0.5 [as suggested by the work of Heavens 
& Drury (1988)] so the small angular size, a < 0.5 
objects from 3C with log 10 P151 ~ 27.5 will not ap- 
pear in our simulations. Particularly in the case of 
the 7C sample, there is a discrepancy between the 
minimum spectral indices observed and those sim- 
ulated and it is unlikely that all of the discrepant 
sources could be explained by invoking synchrotron 
self-absorbtion at these frequencies. In the case of 
the 3C sample, the simulation shown in Figure [Tt] has 
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too high a proportion of sources with linear sizes < 10 
kpc. It is likely that the simplicity of our assumptions 
about the environments of radio sources at these very 
small separations from the AGN are responsible for 
this, in particular the rapidity of source expansion in 
these regimes may be underestimated (as discussed in 
§8.3.2). We do not include in our simulations the ana- 
logues of those sources which we believe to be free-free 
absorbed or synchrotron self-absorbed which may be 
responsible for the absence of these small sources in 
real samples. The bimodality in the simulated spec- 
tral indices is likely to arise in part because we do not 
model any variation in hotspot size in the lifetime of 
an individual source, or in different sources. 

We defer detailed statistical comparison of the dif- 
ferently paramctcriscd simulated datasets and an in- 
vestigation of the constraints on the birth-function in 
the different cosmological models and the purported 
existence of a redshift cutoff in the radio source pop- 
ulation, to a forthcoming paper. 



9. Remarks on interpretation 
9.1. Linear size evolution 

The anti-correlation of linear size with redshift 
arises naturally as a consequence of imposing a survey 
flux-limit on a population of sources each of whose lu- 
minosity falls off with age. For high-power sources at 
high-rcdshift, the P-D tracks are seen to have steep- 
ening slopes throughout the life of the source. This 
arises partly because of the effect of inverse Comp- 
ton losses to the cosmic microwave background (as 
found by Kaiser et al. 1997) and partly because of 
increasing adiabatic expansion losses of successively 
injected particles emerging from a constant pressure 
hotspot into a lobe whose pressure decreases through- 
out the life of the source. In the case of the high- 
est power sources, the injection index itself is steeper 
which causes further steepening of the luminosity of a 
radio source with time. The interception of our light 
cone with a particular object is random throughout its 
lifetime, but high redshift sources in the later stages 
of their lives are increasingly less likely to be detected 
by surveys because of the steep decline of their lumi- 
nosity with age. Thus the only high redshift objects 
we tend to detect are younger (more luminous) ones 
— which being younger, are shorter; hence 'linear size 
evolution' is seen. This means that in a survey, there 
is a strong correlation between youth and redshift in 
the sense that the higher the redshift of a source, the 



earlier in its lifetime we are observing it. Such a de- 
generacy must be taken into account when explaining 
"redshift effects", which we explore in §9.5. 



9.2. Injection indices and the P— a correlation 

The P—a correlation may arise in the way we have 
suggested with the injection spectrum of particles, 
and the steepness of this distribution, governed in 
a way which depends on the jet-power. It is possi- 
ble however, that the acceleration index of immedi- 
ate post-shock particles is not necessarily 0.5 as im- 
plied by first-order Fermi mechanisms, but that this 
spectrum itself is governed by the jet-power in some 
way depending on the detailed physical mechanisms 
of the hotspot. Whatever the details of the underly- 
ing mechanism, an injection spectrum dependent on 
the jet-power is required to give the observed depen- 
dence of spectral index at rest-frame 151 MHz on lu- 
minosity; if the dominant effective break frequency 
were comparable with rest-frame 151 MHz then the 
drastic depletion of the power reservoir of the source 
would have caused it to drop below the flux-limit of 
a typical flux-limited survey. 

Possible evidence that the P-a correlation arises 
in the way modelled in this paper, namely through 
the steepening of the energy distribution of particles 
which are ultimately injected into the lobes by ra- 
diative losses in the enhanced magnetic fields of the 
hotspots of sources with more powerful jets, comes 
from the non-detection of optical synchrotron emis- 
sion associated with the hotspots of very powerful 
radio sources. The number of hotspots which emit 
optical synchrotron radiation is small, for example 
there is a deficit of optical hotspots in very power- 
ful radio sources such as Cygnus A (Roser, Meisen- 
hcimer & Yates 1996), whereas optical hotspots are 
detected in low-power radio sources (Meisenheimer et 
al. 1989). If it is the case, as we have modelled in 
this paper, that more powerful sources have higher 
magnetic fields in their hotspots than do less power- 
ful sources, then the absence of optical synchrotron 
emission from the hotspots of the former arises sim- 
ply because the spectrum of the hotspot emission is 
so steep towards this frequency regime. 
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Fig. 16. — Figure showing real 3C data, for FRIIs and DAs, assuming J7m = 1 & £l\ = 0, for the six planes in 
[PDza] parameter space. Solid circles indicate objects at z > 2, crosses indicate objects with 2 > z > 1.5, open 
circles indicate 1.5 > z > 0.8 and dots indicate sources with z < 0.8. 
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Fig. 17. — Figure showing simulated 3C data, assuming Om = 1 & 17 a = 0, for the six planes in [PDza] parameter 
space and additionally a plot of the jet-power Q a vs rest-frame luminosity at 151 MHz and a plot of the age of each 
source at the point we observe it against redshift. Solid circles indicate objects at z > 2, crosses indicate objects 
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Fig. 19. — Figure showing simulated 7C data, assuming Om = 1 & 17 a = 0, for the six planes in [PDza] parameter 
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If the P—a correlation arises in a way resembling 
that which we have modelled here, then it follows 
that hotspots with more powerful magnetic fields be- 
ing fed by jets which transport higher powers, should 
themselves have steeper spectra than the hotspots in 
sources with lower jet-powers. Such a finding has 
been pointed out by Scheuer (1998), that the pairs of 
hotspots in powerful radio sources have steeper spec- 
tra than those in less powerful sources. 

Previous work [e.g. Baldwin 1982] has shown that 
simple models of radio source evolution predict far 
more high P, large D sources than are actually ob- 
served. We are able to generate a distribution of 
points on the P-D plane which does not give rise to 
such a bulge. This is in part because the higher Q Q 
sources have steeper spectra and so are more easily 
lost below the survey flux- limits; this is in addition 
to the effect whereby the sample selection process it- 
self discriminates against sources which have suffered 
large adiabatic losses. 

We consider the question posed by Jenkins & 
Scheuer (1976) 'what docks the tails of radio source 
components?'. Jenkins & Scheuer found no evidence 
that the tails (i.e. bridge towards the core) of radio 
sources were longer at low-frequency than at high- 
frequency, as might be expected if synchrotron losses 
play a major part in limiting the lengths of tails. We 
contend that a break in the energy spectrum of par- 
ticles injected into the lobe is an increasingly rele- 
vant factor in depleting the bridges of radio sources 
in increasingly powerful sources. An alternative ex- 
planation is that the backflow speeds in these sources 
are too slow for lobe emission to ever reach the core; 
a possible test between these two explanations is by 
consideration of the polarisation patterns in the lobes 
of these sources (Miller 1985). Inverse Compton losses 
arising from the AGN's own radiation field will also 
deplete the lobe emission close into the core (see 



(in) the energy distribution of the particles injected 
into the lobes as discussed in §8.4.2 and (iv) the birth 



§8.1.6) 



9.3. The tightness of the P—z correlation 

Application of a flux-limit to a survey means that 
at the lower luminosity edge there will be a tight cor- 
relation between luminosity and redshift. Just how 
tight the upper edge of the correlation is depends on 
(i) the steepness and/or curvature of the distribution 
in jet-powers (see equation |38| ) (ii) the proximity to 
the upper-limit of the distribution in jet-powers (note 
how the 3C P—z distribution is much tighter than that 
of the 7C distribution comparing Figures HI and [H), 



function itself. The lack of high redshift sources in 3C 
may arise in part because of an intrinsic maximum jet- 
power (= maximum black hole mass) in these objects 
(as discussed in Willott et al. 1998b). 

9.4. Dependence of spectral index on linear 



The tendency for radio sources with larger linear 
sizes to have steeper spectra, as seen in Figure |l is not 
hard-wired into our simulations, yet it does appear 
(e.g. Figure |l9|). We interpret this tendency as be- 
ing due to the decrease in magnetic field as the lobes 
expand. For a fixed observing frequency, the Lorentz 
factor of particles contributing most of their emission 
at this frequency will be higher for a lower magnetic 
field (see equation [|). Thus particles with a higher 
Lorentz factor are required for the chosen emitting 
frequency and given the power-law exponent of typi- 
cal energy distributions, the number of high Lorentz 
factor particles is smaller than the number of less en- 
ergetic particles. Moreover, the adiabatic expansion 
process itself, while preserving the shape of the spec- 
trum, will shift to lower frequencies any features in 
the spectrum such as a break frequency. 

9.5. Ramifications of the 'youth redshift link' 

Each of the Figures [l7] and |l9| show (bottom right 
hand corners) plots of the ages of sources in the given 
simulated sample versus the redshifts. There is a 
very steep (negative) gradient seen in the distribution 
of points on these plots: objects close to a redshift 
of 2 have ages approaching two orders of magnitude 
smaller than objects in the local Universe; the corre- 
spondence in the simulated fainter sample (Figure |l9|) 
is the same but less pronounced. We thus conclude 
that in any study of the importance of redshift to a 
given phenomenon, it is essential to take into account 
the fact that the higher the redshift at which an object 
is observed, the younger (on average) it is. Two such 
issues where this is an important consideration are (i) 
the higher distortion in radio structure of Barthcl and 
Miley's (1988) high-redshift sample of steep-spectrum 
quasars compared with their low-redshift sample and 
(ii) the alignment effect [McCarthy et al. 1987; Cham- 
bers, Miley & van Breugel 1987]. For case (i) Barthel 
& Miley suggested that at higher redshifts the envi- 
ronments of radio sources were more dense and in- 
homogeneous than at low redshift. In case (ii) it has 
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been proposed that the origin of the optical and infra- 
red light which is aligned with the radio axes arises ei- 
ther because of star-formation stimulated by the pas- 
sage of the jet through the medium [McCarthy et al. 
1987; Rees 1989] or by quasar light escaping along the 
radio axes scattered into the line-of-sight by dust or 
electrons [Barthel 1989; Fabian 1989]. 

Sources which are younger may have the inter- 
actions of their jets occuring much closer in to the 
host galaxy where there is a higher density and likely 
greater inhomogeneity in the ambient matter; thus it 
is possible that star-formation is more easily triggered 
here than at distances further out from the central 
engine and possible that the radio structures seen at 
these times are more distorted than those seen much 
later in the life of a radio source. Although a young 
source is not necessarily a short source (consider the 
derivative of equation ^) , in the case of sources which 
are young but quite long on account of their high 
jet-power Q a , the higher jet-power may well be more 
efficient at triggering star formation. 

We note Dey et al.'s (1997) spectropolarimetric 
study of 4C41.17 presents a compelling case for jet- 
triggered star formation in this young, high jet-power, 
radio galaxy at z — 3.8. 

Moreover if it is the case that mergers of galax- 
ies trigger radio galaxy or quasar activity (see e.g. 
Sanders et al. 1988), then observations of young radio 
sources (and thus those at high redshift) will be re- 
vealing the emergence of a radio source through recent 
merger products. Such an environment is more likely 
to be favourable to the triggering of star-formation 
than the rather more 'steady-state' environment of a 
radio source which has been in existence for longer 
than 10 s years. 

Consistent with this picture is the finding of Best, 
Longair and Rottgering (1996) that the smallest sources 
in their sample of z ~ 1 radio galaxies (all with very 
similar luminosities) are those which are most aligned 
with optical emission. Best et al. (1996) remarked 
that the sequence of changing optical aligned struc- 
ture with increasing radio size could be naturally in- 
terpreted by comparing it with different phases of the 
interaction of the radio jets with the interstellar and 
intergalactic medium as the radio sources age. 

We have not needed to invoke any systematic 
change in the environments of radio sources with red- 
shift in order to reproduce good approximations to 
the six planes in [PDza] space from our observations. 



However, we note that Garrington & Conway (1991) 
found that there was a tendency for depolarisation 
to be higher in sources with a higher P (and hence 
z, since their sample mainly comprised objects in the 
3C sample or similar samples) . It is possible that this 
is related to the above discussion: objects in their 
sample with higher P or z which are younger may 
well be in much more recently merged environments 
with the consequence that inhomogeneous regions of 
higher ambient density will more readily depolarise 
polarised intensity radiated from the lobes in these 
environments. 

Note that a result of this work is that we are un- 
avoidably driven to the conclusion that the high red- 
shift sources we find in a radio survey are inevitably 
a younger population than their low redshift counter- 
parts. Comparison of the simulated samples in Fig- 
ures [It] and [l^ suggest that the strong dependence 
of age on redshift is alleviated significantly by studies 
using more than one complete sample e.g. the new 7C 
redshift survey in tandem with the 3C sample. 

Even if some details of our model are not fully cor- 
rect and for example, the inclusion of a source into 
a survey is governed much more by its low-frequency 
emission being dominated by its hotspots than we in- 
fer to be the case on the basis of current data, then in 
order to avoid simulating an excess of very powerful 
large linear size radio sources over that which is seen, 
it would be necessary to invoke that sources at higher 
redshift have their beams switched on for shorter 
times. The "youth-redshift degeneracy" would then 
be truely hard-wired into the model. As it is, in our 
model the "youth-redshift degeneracy" arises because 
older sources do not survive above the survey flux- 
limit because the spectral losses due to (especially) 
inverse Compton and also adiabatic expansion and 
synchrotron losses deplete the luminosity of an ob- 
ject at the relevant emitted frequency [(l+z)x survey 
frequency] . 

The model we have presented here, which repro- 
duces the nature of radio sources which are included 
in low- frequency selected radio surveys, does not need 
to invoke the 'death' of radio sources. This gives us 
a distribution in the ages of the sources which make 
it into our samples but does not require any strong 
statement about their lifetimes. It has been suggested 
that radio sources which are on death row should be 
seen as a population of sources which simply appear 
as widely separated pairs of "naked hotspots". We 
suggest that this is unlikely: the inclusion of such 
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objects into low-frequency (e.g. 151 MHz) surveys is 
most unlikely to be via their hotspots because, as dis- 
cussed in §8.1.5 hotspots do not radiate strongly at 
these low frequencies; although such hotspots will ra- 
diate at frequencies like 1.4 GHz [e.g. as in the FIRST 
survey of Becker et al. (1998)] it is most unlikely that 
emission will be found to link apparently unrelated 
isolated hotspots. 



10. Summary of conclusions 

New observational results 

The population of radio sources which are detected 
by a given survey cannot be naively identified with the 
real population of radio sources out there: our study 
of a number of complete samples of radio sources se- 
lected at frequencies close to 151 MHz, with a cov- 
erage of the luminosity-redshift plane which is sub- 
stantially improved over previous studies, has led us 
to conclude that the principal dependences between 
properties of these radio sources selected in complete 
samples are as follows: 



1. The rest-frame spectral index at low- frequency 
depends on the luminosity of a source [P—a cor- 
relation). 

2. The rest- frame spectral index at high-frequency 
(GHz) depends on the redshift of a source. 

3. Observed- frame spectral indices are also gov- 
erned by the redshift of a source. 

4. Sources with larger physical sizes have steeper 
spectra {D-a correlation). 

5. Sources at higher redshift have smaller physical 
sizes (D-z anti-correlation). 



New views of radio sources 



6. We have presented an intuitive picture for the 
random interception of radio sources with our 
light-cone, which does not necessitate the use 
of assumptions concerning the dependences of 
the lifetimes of radio-sources. 

7. The P-a correlation (result 1, listed above) 
can be understood if the hotspot is viewed as 
the governor of the energy distribution injected 
into the lobes of radio sources, and if the jet- 
power determines the magnetic field density in 
the hotspots. 

8. The D-z anti-correlation (result 4) arises be- 
cause of a more fundamental dependence, namely 
that the high-redshift sources in a given sample 
will be substantially younger than their low red- 
shift counterparts, an effect which we call the 
'youth redshift degeneracy'. 
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Table 1: This table lists the names, redshifts, pro- 
jected linear sizes, flux densities at 178 MHz [on the 



sca le of Ba ars et al. (1977)] and at 38 MHz [ Hales et 



al. |(1995)| and spectral indices between 38 and 178 
MHz of those sources which would not have been se- 
lected were the sample criteria to have been 5 > 60° 
and 5*38 > 25, where S38 is the flux density in Jy at 
38 MHz. 



Name 


z 


Size 
(kpc) 


5*178 

(Jy) 


5*38 

(Jy) 


"38 


3C 268.3 


0.371 


9 


11.7 


12.3 


0.0 


3C343 


0.998 


9 


13.5 


6.3 


-0.5 


3C 343.1 


0.75 


3 


12.5 


9.4 


-0.2 
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Table 2: Table of radio sources from the sample of Allingt on- Smith (1982) selected at 408 MHz in a region of sky 
exactly overlapping with part of the 151-MHz 6C sample (§ |3.2| ) but not common to the 6C sample. [Objects which 
would be too bright at 151 MHz to lie within the flux-limits of the 6C sample are not listed here.] We list their 
names, redshifts, projected linear sizes in kpc, their flux densities at 151 MHz and 408 MHz and their spectral 
indices between 151 MHz and 408 MHz. Those sources below the horizontal rule are those whose flux densities 
at 408 MHz fall below 1.1 Jy; this flux- limit corresponds to those sources with a flux density of 2 Jy at 151 MHz 
and a spectral index of 0.6. The angular sizes of B3 0822+39, B2 0927+35, B3 1056+39, B3 1106+37, B2 1129+37, 
B3 1207+38 and B2 1225+36 were taken from Law-Green et al. (1995); the angular size of B3 1101+38 was taken 
from the NVSS survey (Condon et al. 1998) and all the remainder were taken from the FIRST survey (Becker et 
al. 1998). The redshifts indicated by a | are quoted from Allington-Smith et al. (1988). The redshifts marked by *f 
are quoted from Allington-Smith, Lilly & Longair (1985). The redshifts for sources marked with a single asterisk 
have been measured by S. Rawlings. The redshifts marked with two asterisks are quoted from Eales & Rawlings 
(1996); note that the redshift of B2 1132+37 is highly uncertain. The redshift of B3 1101+38 is quoted from Weiler 
& Johnson (1980). The redshift of B2 1204+34 is quoted from Allington-Smith (1984). The redshift of B2 1225+36 
is quoted from Xu et al. (1994). The redshift of B2 0927+35 was not available to us and so we considered two 
limiting values for this. The 151 MHz flux densities of all these objects were taken from the 6C survey of Hales et 
al. (1988) and Hales et al. (1993) which are complete to ^ 0.2 Jy. For those sources whose names begin with B2, 
their flux densities at 408 MHz are taken from the B2.3 catalogue (Colla et al. 1973) and for those sources whose 
names begin with B3, their flux densities at 408 MHz are taken from the B3 catalogue (Ficarra et al. 1985). 



Name 


z 


Size 


<Sl51 


S408 


a 408 
"151 








(kpc) 


(Jy) 


(Jy) 




B3 0822+39 


1.2* 


8.6 


0.50 


1.76 


-1.3 


B2 0835+37 


0.396* 


32.0 


1.62 


1.26 


0.3 


B2 0927- 


-35 


f 0.100 1 
\ 1.000/ 


r n.i 1 

\ 38.3 J 


. 1.70 


1.15 


0.4 


B3 1049- 


-38 


1.018* 


8.5 


0.99 


1.24 


-0.2 


B31101- 


-38 


0.030 


41.4 


1.70 


1.14 


0.4 


B31106- 


-37 


2.290" 


5.6 


1.70 


1.29 


0.3 


B21129- 


-37 


1.060* 


162.7 


1.30 


1.20 


0.1 


B21132- 


-37 


2.880** 


7.4 


1.10 


1.52 


-0.3 


B31201- 


-39 


0.445** 


6.8 


1.90 


1.11 


0.5 


B3 1056+39 


2.200* 


96.2 


1.70 


0.82 


0.7 


B21204- 


-34 


0.11 


105.0 


1.97 


1.01 


0.7 


B3 1207- 


-38 


0.790** 


221.6 


1.90 


0.94 


0.7 


B21225- 


-36 


1.975 


0.6 


0.50 


1.03 


0.7 


B21245- 


-34 


0.409* 


39.0 


1.80 


1.06 


0.5 



44 



Table 3: The results of the statistical analysis, run for three different sets of cosmological parameters, on all 
members of the complete samples used in the analyses. 



FRIIs and DAs fi M = & fJ A = 





P 


D 


z 


a 


rp a \zD = 0.34 
cr = 5.83 


1~Da\Pz = 

a = 


0.48 

8.59 


r za\PD — 

a = 


-0.32 
-5.34 


z 


rp z \aD = 0.75 
a = 16.01 


rDz\Pa = - 
(T = - 


0.06 
1.01 






D 


rpD\za = -0.13 
a = -2.09 












FRIIs and DAs ft M 


= i & n A = 







P 


D 


z 


a 


rp a \zD = 0.34 
a = 5.75 


fDa\Pz = 

a = 


0.48 
8.58 


r za\PD — 

a = 


-0.24 
-4.05 


z 


rp z \aD = 0.65 
cr = 12.62 


rDz\Pa = - 
cr = - 


0.19 
3.06 






D 


fpD\za = -0.13 

cr = -2.09 












FRIIs and DAs 


1~2m = 


o.i & n A = 


0.9 




P 


D 


z 


a 


rp a \zD = 0.33 
cr = 5.65 


TDa\Pz = 

a = 


0.48 
8.46 


r za \PD = 
a = 


-0.32 
-5.40 


z 


rp z \aD = 0.76 
cr = 16.24 


r Dz\Pa = - 
cr = - 


0.05 
0.85 






D 


rpo\za = -0.12 

cr = -1.88 
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Table 4: The results of the statistical analysis, run for three different sets of cosmological parameters, on just the 
sources with certain FRII structure from the complete samples. 





FRIIs 


n M = o & n 


A = 






P 


D 


z 


a 


rp a \zD = 0.33 
a = 4.91 


TDa\Pz = 0.37 

cr = 5.58 


r za\PD — 

a = 


-0.31 
-4.64 


z 


rp z \ aD = 0.74 
a = 13.63 


rDz\Pa = -0.03 
cr = -0.39 






D 


rpD\za = -0.19 
a = -2.85 










FRIIs 


n M = i & si 


A = 






P 


D 


z 


a 


rp a \zD = 0.32 
a = 4.86 


TDa.\Pz = 0.37 

cr = 5.57 


r za\PD — 

a = 


-0.24 
-3.59 


z 


rp z \aD = 0.61 
cr = 10.32 


rDz\p a = -0.18 
cr = -2.57 






D 


rp D \ za = -0-19 
cr = -2.81 










FRIIs 


r2 M = o.i & o A = 0.9 






P 




z 


a 


rp a \zD = 0.32 
cr = 4.75 


?"Da|P Z = 0.36 

(J = 5.44 


r za \PD = 
a = 


-0.31 
-4.67 


z 


r P z\ aD = 0.74 
cr = 13.87 


TDz\p a = -0-02 
cr = -0.25 






D 


rpD\za = -0.18 
cr = -2.63 
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Table 5: The results of the statistical analysis, run for three different sets of cosmological parameters, on the 3C 
sample alone. 



n M — & Cl A — 





P 


D 


z 




a 


rp a \zD 

a 


= 0.02 
= 0.20 


TDa\Pz = 0.56 

cr = 7.42 


'ZOL 


PD — 
a = 


0.01 
0.12 




z 


rp z \uD 
a = 


= 0.97 
= 25.24 


rDz\Pa = -0.05 
cr = -0.63 








D 


rpD\za = 

a = 


-0.03 
-0.38 














fJ M = 1 & = 








P 


D 


z 




a 


rp a \zD 

a 


= 0.02 
= 0.28 


td&\Pz = 0.56 
cr = 7.42 


'ZOL 


PD — 

a = 


0.03 
0.33 




z 


rp z \aD 

a = 


= 0.96 
= 22.77 


Tdz\Po. = -0-12 

cr = -1.39 








D 


r PD\za — 

a = 


-0.02 
-0.20 














J1m = 


= o.i & n A = 0.9 








P 


D 


z 


a 


rp a \zD 

a 


= 0.03 
= 0.35 


tdu\Pz = 0.55 
o- = 7.32 


^za 


PD — 

a = 


-0.01 
-0.13 


z 


r Pz\aD 

a = 


= 0.98 
= 26.20 


rDz\p a = -0-03 
cr = -0.33 








D 


TpD\za = 

a = 


-0.04 
-0.47 
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Table 6: The parameters n and m from equation [l] which have been fitted to the data for the three different 
assumed cosmologies. 

n M = o & 12 A = o 

FRIIs and DAs 



n[rdz\pa] = 0.40 ±0.30 

n[rdz\p] = 1.52 ±0.30 

n[rdz\a] = 0.86 ± 0.20 

m[rpd\za] = 0.02 ±0.15 

m[rpd\z] = 0.01 ±0.20 

m[rpd\a] = 0.07 ±0.20 

n M = i & 12 A = o 

FRIIs and DAs 



n[rdz\pa] = 0.96 ± 0.25 

n[rdz\p] = 1.87 ±0.30 

n[rdz\a] = 1.31 ±0.20 

m[rpd\za] = 0.01 ±0.15 

m[rpd\z] = 0.09 ±0.30 

m[rpd\a] = 0.26 ± 0.20 

n M = o.i & n A = 0.9 

FRIIs and DAs 



n[rdz\pa] = 0.35 ± 0.30 

n[rdz\p] = 1.41 ±0.35 

n[rdz\a] = 0.76 ± 0.20 

m[rpd\za] = 0.02 ± 0.20 

m[rpd\z] = 0.07 ±0.30 

m[rpd\a] = 0.18 ±0.20 
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